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Abstract. The ternary Goldbach conjecture states that every odd number 
n > 7 is the sum of three primes. The estimation of sums of the form 
<x e(ap), a = a/q + 0(l/q 2 ), has been a central part of the main ap- 
proach to the conjecture since (Vinogradov, 1937). Previous work required q 
or x to be too large to make a proof of the conjecture for all n feasible. 
, The present paper gives new bounds on minor arcs and the tails of major 

>>: 



arcs. For q > 4 • 10 6 , these bounds are of the strength needed to solve the 
ternary Goldbach conjecture. Only the range q € [10 5 ,4 • 10 6 ] remains to be 
checked, possibly by brute force, before the conjecture is proven for all n. 

The new bounds are due to several qualitative improvements. In particu- 
lar, this paper presents a general method for reducing the cost of Vaughan's 
■ identity, as well as a way to exploit the tails of minor arcs in the context of 

the large sieve. 
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1. Introduction 



The ternary Goldbach conjecture (or three-prime conjecture) states that every 
odd number > 7 is the sum of three primes. I. M. Vinogradov [Vin37] showed 
in 1937 that every odd integer larger than a very large constant C is indeed the 
sum of three primes. His work was based on the study of exponential sums 



and their use within the circle method. 

Unfortunately, further work has so far reduced C only to e 3100 ([LW02]; see 
also [CW89] ). which is still much too large for all odd integers up to C to be 
checked numerically. The main problem has been that existing bounds for (jl.ip 
in the minor arc regime - namely, a = a/q + 0(l/q 2 ), gcd(a, q) = 1, q relatively 
large - have not been strong enough. 

The present paper gives new bounds on smoothed exponential sums 

(1.1) S v (x, a) = A(n)e(an)r)(n/ x) . 




These bounds are clearly stronger than those on smoothed or unsmoothed expo- 
nential sums in the previous literature, including the bounds of |Tao] . (See also 
work by Ramare [RamlO] .) 

In particular, on all arcs around a/q, q > 4- 10 6 , the bounds are of the strength 
required for a full solution to the three-prime conjecture. The same holds on the 
tails of arcs around a/q for smaller q. 

Taken together with existing computational work [Plall] . this means that all 
that is needed to complete the proof of the three-prime conjecture is a bound for 
a finite number of arcs, namely, those around a/q, 10 5 < q < 4 • 10 6 . Moreover, 
these arcs are made rather narrow by the estimates given here on tails of arcs. 
Giving a bound for this narrow range of g's is a task that could be done either 
by a finite amount of computation - namely, an extension of the computations in 
[Plall| to all q < 4 • 10 6 - or by a further improvement of bounds on exponential 
sums. 

The quality of the results here is due to several new ideas of general applica- 
bility. In particular, §4.1l introduces a way to obtain cancellation from Vaughan's 
identity. Vaughan's identity is a two-log gambit, in that it introduces two con- 
volutions (each of them at a cost of log) and offers a great deal of flexibility in 
compensation. One of the ideas in the present paper is that at least one of two 
logs can be successfully recovered after having been given away in the first stage 
of the proof. This reduces the cost of the use of this basic identity in this and, 
presumably, many other problems. 

We will also see how to exploit being on the tail of a major arc, whether in the 
large sieve (Lemma 14.31 Prop. 14. 5p or in other contexts. 

There are also several technical improvements that make a qualitative differ- 
ence; see the discussions at the beginning of $3] and $H Considering smoothed 
sums - now a common idea - also helps. (In the general context of the circle 
method, smooth sums go back to Hardy-Littlewood; in recent work on the current 
problem, they appear in [Tao].) 
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1.1. Results. The main bound we are about to see is essentially proportional to 
((log q) I W 4>(q)) • x. The terms <5o, *i below are there to improve the bound when 
we are on the tail of an arc. 

Main Theorem. Let x > x$, xq = 3.1 • 10 28 . Let S^(a, x) be as in with r/ 

defined in Let a = a/q + 5/x, q < Q, gcd(a, q) = 1, \5/x\ < 1/qQ, where 

Q = x °- 7 /8. Lfq< 2x - 3 , then 

._, ..^ 3.2423 + R xAg log 28 x q 1.6943 + 0.1108 log 2S iq 
\b v {a, x)\ < • x H • x 

( 12 ) f min log x + \ log q + 0.5683, f log q + 17. 



S&(q) 



x 



+ 0.0045— T + 3.336x 



0.85 



where 
(1.3) 



5 = max(l, |tf|/4), *l = if rl/ ^j 5 j- 4 * 2 =max(l 



7T 



i^. t = 0.4091 + 0.28151 log 1 + 



\ 1*1/2 i/i*i>4, ""-"T'vr 

log 50^ \ 



> 2x - 3 , i/ien 



„0.85 



|5^(a,x)| < 0.2974(loglogx)(log3;)x t 

For x > xq and g < 2x 03 , R x>q < R xo2x o.3 = 0.96188... and R x ,s iq < 
^xo,4xO-3 = 1.11890. . . . If q ~ 4 • 10 6 and |*| < 4 (likely to be the worst sit- 
uation in practice, given major-arc bounds) then R Xj s iq < R XOy 2-io 6 = 0.6238 .... 

The classical choice for r\ in (jl.ip is ry(t) = 1 for t < 1, r/(i) = for t > 1, 
which, of course, is not smooth, or even continuous. We use 

(1.4) r](t) = % (t) = 4max(log2 - | log 2*|,0), 

as in Tao [Tao], in part for purposes of comparison. (This is the multiplicative 
convolution of the characteristic function of an interval with itself.) Nearly all 
work should be applicable to any other sufficiently smooth function rj of fast 
decay. It is important that rj decay at least quadratically. 

1.2. History. The following notes are here to provide some background; no claim 
to completeness is made. 

Vinogradov's proof [Vin37| was based on his novel estimates for exponential 
sums over primes. Most work on the problem since then, including essentially 
all work with explicit constants, has been based on estimates for exponential 
sums; there are some elegant proofs based on cancellation in other kinds of sums 
([HB85], [IKQH §19]), but they have not been made to yield practical estimates. 

The earliest explicit result is that of Vinogradov's student Borodzin (1939). 
Vaughan |Vau77] greatly simplified the proof by introducing what is now called 
Vaughan's identity. 

The current record is that of Liu and Wang [LW02J: the best previous result 
was that of |CW89| . Other recent work falls into the following categories. 

Conditional results. The ternary Goldbach conjecture has been proven under 
the assumption of the generalized Riemann hypothesis [DEtRZ97j. 
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Qo 


\S v (a/q,x)\ , 

x nere 


\S v (a/q,x)\ ry— l 


1U 


U.UoOoO 


U. 40000 


10 6 


0.03165 


0.17515 


2 • 10 6 


0.02350 


0.13209 


3 • 10 6 


0.01977 


0.11220 


4- 10 6 


0.01751 


0.09999 


10 7 


0.01194 


0.06934 



Table 1. Upper bounds here and in |Taoj on x 1 \S r) (a/q, x)\ for 
q > qo, \S\ < 4, x = 3.1 ■ 10 28 . The trivial bound is ~ 1. 



Ineffective results. An example is the bound given by Buttkewitz [ButllJ. The 
issue is usually a reliance on the Siegel-Walfisz theorem. In general, to obtain 
effective bounds with good constants, it is best to avoid analytic results on L- 
functions with large conductor. (The present paper implicitly uses known results 
on the Riemann £ function, but uses nothing at all about other L-functions.) 

Results based on Vaughan's identity. Vaughan's identity jVau77] greatly sim- 
plified matters; most textbook treatments are by now based on it. The minor-arc 
treatment in [Taoj updates this approach to current technical standards (smooth- 
ing), while taking advantage of its flexibility (letting variable ranges depend on 
I)- 

Results based on log-free identities. Using Vaughan's identity implies losing a 
factor of (logx) 2 (or (logg) 2 , as in |Tao] ) in the first step. It thus makes sense to 
consider other identities that do not involve such a loss. Most approaches before 
Vaughan's identity involved larger losses, but already [Vin371 §9] is relatively 
economical, at least for very large x. The work of Daboussi [Dab96j and Daboussi 
and Rivat [DROl] explores other identities. (A reading of [DROl] gave part of the 
initial inspiration for the present work.) Ramare's work [RamlOj - asymptotically 
the best to date - is based on the Diamond-Steinig inequality (for k large). 

* * * 

The author's work on the subject leading to the present paper was at first 
based on the (log- free) Bombieri-Selberg identity (k = 3), but has now been 
redone with Vaughan's identity in its foundations. This is feasible thanks to the 
factor of log regained in E j4.ll 

1.3. Comparison to earlier work. Table Q] compares the bounds for the ratio 
\S v (a/q,x)\/x given by this paper and by |Tao] for x = 3.1 • 10 26 and q > qo, 
where qo is 10 5 , 10 6 or 4 • 10 6 . We are comparing worst cases: 4>(q) as small as 
possible (q divisible by 2-3-5 • • • ) in the result here, and q divisible by 4 (implying 
4a ~ a/(g/4)) in Tao's result. The main term in the result in this paper improves 
slowly with increasing x; the results in [Taoj worsen slowly with increasing x. 

We are focusing on a comparison with [Tao] for two reasons: (a) the results 
in [DROl] are unfortunately worse than the trivial bound in this range, and 
(b) Ramare's result [Rami 01 Thm. 3] is not applicable in this range, since its 
condition logg < (l/50)(log x) 1 / 3 is not satisfied. 

The qualitative gain with respect to [Taoj is about log(g) \/4>{q)/q, which is 
~ log(g) / \J e~< (log log q) in the worst case. 

As we shall see, Table Q] suggests that the main theorem in the present pa- 
per will be enough to prove the ternary Goldbach conjecture once the zeroes of 
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Dirichlet L- functions of conductor 10 5 < q < 4 • 10 6 are checked up to a moderate 
height. (Since the bound in the main theorem depends on 1/v/ (f>(q) rather than 
1/sqrtq, it would actually be enough to check all conductors 10 5 < q < 10 6 and 
some conductors 10 6 < q < 4 • 10 6 , e.g., those that are divisible by 30.) 

1.4. Existing numerical work. The numerical inputs we can count on are as 
follows 

• every odd number n < tlq is the sum of three primes, where 

n = 6.15697 • 10 28 . 

This is a moderately cautious choice for no, in that there are statements 
in the literature that would enable one to take no higher; see Appendix 

• L-functions of conductor q < 10 5 have no zeroes outside the critical line 
up to height \Q(z)\ < W s /q (D. Piatt's thesis [Plallj ). 

The estimates in this paper for large 5, together with D. Piatt's results on 
zeroes of L-functions of small conductor, should be enough to give major-arc 
bounds for q < 10 5 . 

Running Piatt's program to extend his results to all conductors q < Qo up to 
height |3f(^)| < QoT/q with T = 0(Qq) would take time 0(Q 3 log 2 Q); consider- 
ing a single q takes time 0(QT log q log QT) (D. Piatt, personal communication). 
Since the resources needed for Qo = 10 5 were relatively modest, this would seem 
to put the task of checking all q up to 10 6 and some q up to 4 • 10 6 just about on 
the limit of what is nowadays feasible by a purely computational approach. 

1.5. Major and minor arcs: a summary assessment. In the circle method, 
the number of representations of a number N as the sum of three primes is 
represented as an integral over the "circle" R/Z, which is partitioned into major 
arcs 931 and minor arcs m: 

^2 A(ni)A(n2)A(n 3 ) = / (S(a,x)fe{-Na)da 

(15) ni +n 2 +n 3 =N 

(S(a,x)fe(-Na)da + / (S(a, x)fe{-Na)da. 



/an Jm 
The major arcs 9Jt consist of intervals around the rationals a/q, q < Qq. In 
previous worlo, Qo grew with x; in our setup, Qo is a constant. Smoothing 
changes the left side of ([1.5]) into a weighted sum, but, since we aim at an existence 
result rather than at an asymptotic for the number of representations p\ +P2 +P3 
of N, this is obviously acceptable. 

1.5.1. Minor arcs. The main theorem bounds max a6m \S(x, a)\. The ternary 
Goldbach conjecture will be proven once we obtain an upper bound for the minor- 
arc integral If \S(x, a)| 3 da| that is smaller than a lower bound for the major-arc 
integral f m S(x, a) 3 e(— Na)da. 

(See $6]for a more detailed discussion.) 

Traditionally, the upper bound for the minor-arc integral | f \S(x, a)| 3 do:| took 
the form 



(1.6) 



\S(x, a)\ 3 da 



< (max|5(a,x)|oo) • \S\l- 



' ttgm 



^Ramare's work |RamlO| is in principle strong enough to allow Qo to be an unspecified large 
constant. Tao's work [Tap] reaches this standard only for x of moderate size. 
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One can of course improve on this by taking out the major-arc contribution: 



(1.7) 



\S(x, a)\ 3 da 



< (max \S(a, x)\ oa ) ■ ( — / \S(x, ct)\ A da 



an 



More crucially, an idea of Heath-Brown's (based on Montgomery's inequality 
[Mon68| . |Hux72j : see [Taol Lemma 4.6]) allows us, in essence, to multiply the 
term \S\ 2 in (JL6]) or (jl.7jl by a factor that, while only slightly less than 1 for 
x ~ xo = 3.1 • 10 28 and Qo ~ 4 • 10 6 , grows like 0(1/ log x) as x grows. It is this 
that allows us to let Qo be a constant. 

A brief numerical check shows that, together, 

• the factor coming from Heath-Brown, 

• the fact that R x ,& iq in fll.2j) decreases as x increases, 

are enough to make the situation for x > xq, q ~ 4 • 10 6 better when x ^> xq 
than when x ~ xq. Thus we can focus on x ~ xq as the worst case. 

In summary, because x ~ xq is the worst case and because we can reduce a 
little the factor of \S\?. ~ logxo in (|1.6p . a bound of slightly less than 0.02x on 
the minor arcs is necessary and sufficient for the minor-arc contribution to be 
acceptable. By Table [H this corresponds to Qo = 3 • 10 6 ; we choose Qo = A - 10 6 
to allow for losses in the major-arc part of the argument. 

1.5.2. Major arcs. Since Qo is a constant, the treatment of the major arcs is a 
finite problem that can be attacked computationally. The actual estimation of 
major arcs is a subject for a separate paper; we will simply discuss very briefly 
what the actual size of the term being estimated should be. 

We can choose a different smoothing rji for the sum whose £2 norm is taken 
(in (jl.6p or (|1.7j) ). as opposed to the smoothing 770 in the sum whose ioo norm we 
bound in the main theorem. 

For at least moderately large Qo ( as is the case here), the major-arc contribu- 
tion is close to what is believed to be the actual weighted number of representa- 
tions of N: 
(1.8) 

(S ri0 (a,x))(S ril (a,x)) 2 da ~ ^ A(ni)A(n 2 )A(n 3 )r? (n/x)r/i(n/x)?7i(n/x) 

1 



9H 



(p-i) 2 

p>3 v yF ' 



= 1.3203... • (r/o* r/i*??!)^), 

where f3 = N/x. 

Recall that the minor-arc estimate is proportional to |<SL ||, and hence to 1 771 1 2 . 
Hence, we should choose /3 (within a given range) so as to make 

9 ) Vo *Vi *Vi(P) 

Ml 

as large as possible; this should also play a role in our choice of 771. It is easy to 
see that (jl.9p is at most 1 771 * 7/iIqo • | t7q 1 1 / 1 t?x 1 2 < |t?o|i = 1; the optimal value of 
1 is approached when 771 is roughly symmetric (additively, around some positive 
value) and the support of 771 is somewhat broader than that of 770. 

There are other matters in play: it is also desirable to choose an 771 whose 
Mellin transform decays very rapidly, as this minimizes the error term in the 
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major-arcs estimate. This means we should be content with a good enough 771, 
rather than search for an rji that is optimal according to a single criterion. 

1.5.3. Conclusion. The minor-arc estimate, without the substraction of a major- 
arc term as in (jl.7p . is 

15. I 2 

(1.10) 1.163699 • ~ 1-163 . . . • \m\l 

log x 

for Q = 4 • 10 6 and x > x = 3.1 • 10 28 . Subtracting a major-arc term would 
reduce this. 

The major-arc integral is close to 

(1.11) 1.3203... *m * Vl) W). 

If the total losses from an error term here and the ratio (|1.9p are less than 10 
percent, the three-prime theorem will be proven. Larger errors can obviously be 
allowed if Qq is set higher. 

1.5.4. Final remarks. Larger losses will be allowable if q is set larger than 4 • 10 6 
or the effect of the substraction of a major-arc term from (|1.10|) is substantial. 
It is likely possible to omit all odd values of q between 2 ■ 10 6 and 4 • 10 6 from 
consideration, as the terms 4>{q) in (11.21) are then closer to q. A technical idea in 
|Taoj (using 4a instead of a, and summing over odd numbers in some contexts) 
is also likely to help for q odd, or q = 2 mod 4. 
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crucial help and feedback, and to D. Piatt for his prompt and helpful answers. 
He is also much indebted to A. Booker, B. Green, H. Kadiri, T. Tao and M. 
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are also due to B. Bukh, A. Granville and P. Sarnak for their valuable advice. 

Travel and texts were funded in part by the Philip Leverhulme prize. The 
author's work on the problem started at the Universite de Montreal (CRM) in 
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The present work would most likely not have been possible without free and 
publicly available software: Gnuplot, PARI, QEPCAD, SAGE, and, of course, 
LT^X, Emacs, the gcc compiler and GNU/Linux in general. 

2. Preliminaries 

2.1. Notation. Given positive integers m, n, we say m|n°° if every prime divid- 
ing m also divides n. We say a positive integer n is square-full if, for every prime 
p dividing n, the square p 2 also divides n. (In particular, 1 is square-full.) We 
say n is square-free if p 2 \ n for every prime p. For p prime, n a non-zero integer, 
we define v p (n) to be the largest non- negative integer a such that p a \n. 

When we write we mean X^Li> unless the contrary is stated. As usual, /z, 
A, r and a denote the Moebius function, the von Mangoldt function, the divisor 
function and the sum-of-divisors function, respectively. 

As is customary, we write e(x) for e 2mx . The Fourier transform on M is nor- 
malized here as follows: 

/oo 
e{-xt)f(x)dx. 
-oo 

We write \f\ r for the L r norm of a function /. 



H. A. HELFGOTT 



We write 0< (R) to mean a quantity at most R in absolute value. We quote the 
following standard bound from |Tao[ Lemma 3.1]: for a £ M/Z and / : K — > C 
compactly supported, continuous and piecewise C 1 , 



(2.1) 



<min fl/K + Vlx, ,4^ 
\ 2 sm to 



When we need to estimate ^ n /(n) precisely, we will use the Poisson summa- 
tion formula: 

n n 

We will not have to worry about convergence here, since we will apply the Pois- 
son summation formula only to compactly supported functions / whose Fourier 
transforms decay at least quadratically. 

If / is compactly supported and C 2 , then the Fourier transform of / does decay 
at least quadratically: 



/(*) 



(2.2) 



f(x)e(—tx)dx 
tx) 



m e -^dx 



f"(x) 



dx = 0< 



1 



(2-irty 



-2nit 
■\f"h 



' (2vut) 2 

This can still be made sense of when /" is undefined at a finite number of points, 
provided / is understood as a distribution (and /' has finite total variation). Of 
course, \f(t)\ < \f\i, as always. 

For the smoothing function t]q in (|1.4jl . 

(2-3) Mi = l, K|i = 81og2, |77o|l = 48, 

as per |Taol (5.9)-(5.13)]. Similarly, for rji, p (t) = log(pt)r/o(t), where p > 4, 

\m, P \i < log(p)|77o|i = log(p) 

Wiji = 2mA 1 / 2 ) = 21og(p/2) 7? o(l/2) < (8 log 2) log p, 
U " > \ViJi =4Iog(p/4) + |21ogp- 41og(p/4)| + |41og2 - 41o g/ o| 
+ | log p — 4 log 2| + | log p\ < 48 log p. 

Write log + x for max(logx, 0). 

2.2. Bounds on sums of p(m) and A(n). We will need explicit bounds on 
J2n<N l I i. n )/ n an d related sums involving p. The situation here is less well- 
developed than for sums involving A. The main reason is that the complex- 
analytic approach to estimating ^2 n<N p(n) would involve 1/C(s) rather than 
C'(s)/C(s), and thus strong explicit bounds on the residues of 1/C(s) would be 
needed. 

Fortunately all we need is a saving of (logn) or (logn) 2 on the trivial bound. 
This is provided by the following. 

(1) ( Granville- Ramare [GR96], Lemma 10.2) 



(2.5) 



E 

n<x:gcd(n,<j)=l 



p(n) 



n 



< 1 



for all x, q > 1, 
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(2) (Ramare |Ramc| : cf. El Marraki |EM95j . |EM96| ) 



(2.6) 



£ 

n<x 



H(n) 



n 



< 



0.03 
log x 



for x > 1 1815. 
(3) (Ramare [Ramb] ) 

(2.7) £ ^ = O s ( m in( 2 . 

n<x:gcd(n,q)=l 

for all x and all q < x; 

/i(n 



02 



<j){q) \ogx/q 



■1 



x 



n 



log - = 0< 1.4 + 



n 



(2-8) £ 

n<x:gcd(n,g)=l 

for all x and all g. 

Improvements on these bounds would lead to improvements on type I estimates, 
but not in what are the worst terms overall at this point. 

By numerical work carried out by Olivier Ramare and reproduced indepen- 
dently by the author, 



(2.9) 



E 

n<x 



H(n) 



n 



x 



for all real x > 10 10 . (In fact, the stronger bound < l/2y^ holds for all 3 < x < 
7727068587, but not for x = 7727068588 - e. The author's verification of (|23j) 
goes up to 10 11 .) 

We will make reference to various bounds on A(n) in the literature. The 
following recent bound ([Rama], supplemented by a numerical check for very 
small x) will be useful: 



(2.10) 



E 

n<x 



A(n) 



n 



< logx. 



We also use the following older bounds: 

(1) By the second table in |RR96} p. 423], supplemented by a computation 
for 2 • 10 6 < V < 4- 10 6 , 

(2.11) £ A (n) < 1.0004?/ 

n<y 



for y > 2 ■ 10 6 



(2 

(2.12) 



J2 A(n) < 1.03883y 

n<y 

for every y > [ RS621 Thm. 12]. 
For all y > 663, 



(2.13) 



£ A(n)n < 1.03884y, 
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where we use (|2.12p and partial summation for y > 200000, and a computation 
for 663 < y < 200000. It is also true that 

(2.14) ( l °SP) 2 <lv(logy) 

y/2<p<y 

for y > 117: this holds for y > 2 ■ 758699 by [RS751 Cor. 2] (applied to x = y, 
x = y/2 and x = 2y/3) and for 117 < y < 2 ■ 758699 by direct computation. 

2.3. Basic setup. We begin by applying Vaughan's identity |Vau77| : for any 
function i] : M — > M, any completely multiplicative function / : Z + — > C and any 

x > 0, U, V > 0, 

(2.15) ^2 Hn)f(n)e(an)r](n/x) = S/,1 - S I>2 + S n + ^0,00, 

n 

where 

S/,1 = ^2 M( m )/( m ) y~^(l°g n)e(amn)f(n)r](mn/x), 
m<U n 

Si ; 2 = A(d)/(d) j u ( m )/( m ) e(admn)f(n)r](dmn/x), 



(2.16) 



S U = ~£ f(m) 



m>U 



/i(d) A(n)e(amn)f(n)rj(mn/x), 

d>U n>V 



\d\m / 
So,oo = ^2 H n )e(an)f(n)rj(n/x). 



n<V 



The proof is essentially an application of the Mbbius inversion formula; see, e.g. 
[IK04J §13.4]. In practice, we will use the function 



/(") 



j 1 if gcd(n, it?) = 1, 
I otherwise, 

where w is a small positive, square-free integer. Then 

S v (x, a) = S^i — + Sn + So,oo + Sq iW , 
where S v (x,a) is as in (jl.ip and 

So,w = A(ri)e(an)r](n/x). 

n\w 

The sums Si t i, Sip are called "of type I" (or linear), the sum Sn is called "of 
type II" (or bilinear). The sum So is in general negligible; for our later choice of 
V and 77, it will be in fact 0. The sum So >w will be negligible as well. 

3. Type I 

There are here three main improvements in comparison to standard treatments: 

(1) The terms with m divisible by q get taken out and treated separately 
by analytic means. This all but eliminates what would otherwise be the 
main term. 

(2) For large m, the other terms get handled so as to replace a factor of log q 
by a factor of log min(W/U, q). 
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(3) The "error" term 5/x = a — a/q is used to our advantage. This happens 
both through the Poisson summation formula and through the use of two 
successive approximations. 

3.1. Trigonometric sums. The following is a Vinogradov-type lemma (occu- 
pying roughly the middle ground between Lemma 8a and Lemma 8b in [Vin04, 
Ch. I]). 

Lemma 3.1. Let a = a/q + [3/qQ, gcd(a,q) = 1, \/3\ < 1, Q > q. Then, for 
U>1, 

(3.1) V min(V, 1 ^ ^ ] <min( 2U + — log2q,3U + — log2eLM . 

^-^ V smfTran) / V 7r tt / 

y<n<y+q v 7 

In contrast to the standard procedure, we will sometimes choose to bound 
some terms by U and other terms by 1/| sin(7ran)|. This simple change will allow 
us - in essence - to replace a factor of log q by a factor of log U whenever this 
is to our advantage (see (13. 1|) ). We otherwise follow Daboussi and Rivat [DR01, 
Lemma 1], with minor improvements. 

Proof. We write tuq = [y] + [(g + l)/2], j = n — rriQ, so that j ranges in the interval 
{-q/2,q/2\. Then 

a 3 Pj a 3 c- r / .\ 

an = am H h — — = \-Oq + qAj) mod 1, 

q qQ q 

where 5q = am® is fixed and 5\{j) = [3j/qQ satisfies |5i(j)| < 1/2Q < l/2q. 
Write <5o in the form c/q + 62, where c is an integer and l^l < 1/2(7- As j runs 
over the interval —q/2 < j < q/2, the variable r = aj + c modg occupies each 
residue class mod p exactly once. We may assume 82 > without loss of generality 
(since we can replace a by —a if needed). We bound the terms corresponding to 
r = and r = — 1 by U each; we will now consider the other terms. 

One option (advantageous for U large, but always valid) is to bound every term 
min(C7, 1/| sin(</>an)|) by 1/| sin(</>an)| from now on. The terms corresponding to 
r = —k and r = k — 1 (2 < k < q/2) contribute at most 

1 111 

+ - ; ; r < 7 + ■ 



sin \{-k + 1/2 + ^2) sin 3/2 + qS 2 ) ~ sin | (k -1/2) sin ^(k - 3/2) ' 

since x 1— > 1/ sinx is an even function that is convex- up on (0, 00). Hence the 
terms with r 7^ 0,-1 contribute in total 

1^1 1 f q ' 2 1 , 

— 7—, — r + 2 > ^7 t-t < — 7—, — r + 2 / ^—dx, 

sm(7r/2q) 2 ^ /2 sin|(r-l/2) " sin(vr/2g) J 1 sin^x 

where we use again the convexity of x 1— > 1/sinx. Now 

9/2 1 , g, vr g , 2q 

ax < — log cot — < — log ■ 



sin —x tt 2q ir tt 

and since we may assume q > 3 (the statement of the lemma is trivial otherwise) 
we have l/sin(7r/2g) + {2q/n) log(2g/vr) < (2q/ir) log 2q. 

The alternative is to use the bound < U for some terms with r 7^ 0,-1. We 
consider this only when U < 1/ sin(7r/2q). Assume C/ > 1 (or else the estimate 
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< Uq holds and is optimal for trivial reasons). The terms with r 7^ 0,-1 then 
sum up to at most 



(3.2) 



U+ V mm [ U, . n , -. 

^ \ sm- r'-] 

2<r'<g/2 V 9 V 

+ V min f U, . n . - 1 
^— ' \ sm - r — 1 

2<r'<«/2 \ 9 V 



/2 + g5 2 ) 



We choose to bound a term 1/ sin((7r/g , )(r / — 1/2 ± g^)) by {/ if and only if 
1/ sin((7r/g)(r' — lig^)) > EP The number of such terms is < (g/vr) arcsin (1/ZL7")=F 
g52, i.e., (2q/n) arcsin (1/U) in total. (Note we are using (g/vr) arcsin(l/E/) > 
1/2.) Each other term gets bounded by the integral of 1 / sva.{<pa / q) from r' — lit 
q&2 (> (<7/tt) arcsin(l/E/)) to r' ± g#2, by convexity. Thus (|3.2|) is at most 

2a 1 Z" 9 / 2 1 

E/ + —U arcsin — + 2 / ^-^dt 



VT U Jr ar csm(l/l/) sin f 

< E7 H — -E7 arcsin — + 2— log cot f - arcsin — 

7T (7 7T \2 U 

< C/ + ^(l + log(2C/)). 

7T 

A quick comparison shows that U + (2g/7r)(l + log(2E/)) < (2g/7r) log 2g when 
U < 0.2488... times q. Note U < q/A implies U < 1/ sm(ir/2q), so the bound 
U + (2q/n)(l + log(2E/)) is certainly valid for U < q/A. □ 

Lemma 3.2. Let a = a/q + /3/qQ, gcd(a, q) = 1. |/3| < 1, Q > q. Let 1/2 > Hi > 
0. If 1/2 — Hi < q and 1/2 < Q/2, then, for any U > 1, 

V minf U,-. t~ rr ) <minf — log2eg,2E/ + — log2eE7 J 

^— ' \ | sin(7rem) | / \ir it J 



to o\ y^< n <y2 

(3-3) g| n 



< — log 16E7. 

7T 



This is a variant of Lemma [3.2l with the term for n divisible by q set aside. The 
bound (2q/ir) log 2eq is smaller than 2U + {2q/ir) log 2eE7 for {7 of size roughly 
g/3 or larger. 

Proof. For n < 2/2 > a™ equals an/q + f3n/qQ = an/q + 0<(l/2q). We proceed as 
in the proof of Lemma 13-H only without excluding the terms labelled r = 0, —1 
there; the sum we obtain is at most 

(3.4) 2 min ( E/". 



r=l \ <? 



s i n ZL( r _ 1/2) 



The estimate without EE's is 2/sin(7r/2(7) + (2q/ir) log(2g/7r) < (2q/ir) log 2eg. 
For the estimate with E/'s, proceed as before, using the convexity of 1/sin and 
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choosing U whenever r < 1 + (q/n) arcsin(l/?7): (|3.4p is at most 

rq/2 1 



2u(l + - arcsin(l/C/) N ) + 2 f 



( g /7r)arcsin(l/l/) sm f * 



< 2C7 H — - ( U arcsin h log cot — arcsin — 

7T y [7 2 U 

<2U+ — log 2eC/. 

7T 

To obtain the last line of (J373J), note that, if q > (16/2e)J7, then 2f7 < y§ < 
(2/vr)(logl6-log2e). □ 

3.2. Type I estimates. Our main type I estimate is the following. One of the 
main innovations is the manner in which the "main term" (m divisible by q) is 
separated; we are able to keep error terms small thanks to the particular way in 
which we switch between two different approximations. 

Lemma 3.3. Let a = a/q + 5/x, gcd(a, q) = 1, \S/x\ < 1/q 2 , q < x. Let rj be 
continuous, piecewise C 2 and compactly supported, with rj" G L\ . Let w G Z + be 
square-free. Then, for D < x/w and any e £ (0, 1] such that e/(l + e) > 2\5q\/x, 



m<D n 
gcd(m,w>)=l gcd(n,w)=l 



/i(m) e(amn)t](mn/x) 

n 

gcd(n,u>)=l 

= x(J,(q)(/)(w/q + ) _ //(to) 

q w/q + m 

' m<M/g- 

gcd(m,wq-)=l 

+ { xq/w "U"^J k 11 

, , . / q , Dw ( , . , 2ecx , 16ecx \ 
+ r?i • ^log2eg + 2(1 + e) log — - + log — — 

\7T 7T \ D Dw J 

, ,,x / 3wcx 3(l + e)x, Dw\ 

+ (2|ml + WI)(-g- + -L_J_ log _j, 

w /iere c = 2|7?|i/|r/'|i + 1, Q = |z/l<5<?IJ, ^ = min(LDj , Q/2w), q + = gcd(q,w) 
and q~ = q/q + ■ 

A few words on error terms and "recommended usage". The last two lines in 
(|3.5p are usually very small. The term D log((4ecx) / D) is independent of our 
choice of approximation; it forces us to have D smaller than x times a constant 
if we want an error smaller than x. If a, q are chosen so that \a — a/q\ < 1/qQo 
for some Qo, then Q = [x/\5q\\ > Qq. If 5 = 0, we set Q = oo. 

Proof. Since Q = \x/\8q\\, we have a = a/q + 0<(l/qQ) and q < Q. (If 5 = 0, we 
omit what follows.) Let Q' = [~(l + e)Q] > Q + l. Then a is not a/q + 0<(l/qQ'), 
and so there must be a different approximation a'/q', gcd(a',q') = 1, q' < Q' 
such that a = a'/q' + 0<(l/q'Q') (since such an approximation always exists). 
Obviously, \a/q — a! /q'\ > 1/qq', yet, at the same time, \a/q — a'/q'\ < 1/qQ + 
1/q'Q' < 1/qQ + 1/((1 + e)q'Q). Hence q'/Q + q/{{\ + e)Q) > 1, and so q' > 
Q-q/(l + e) > (e/(l + e))Q. (Note also that (e/(l + e))Q > (2|tfg|/a;)- [a;/5gj > 1, 
and so </ > 2.) 
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By Mobius inversion, 



fj,(m) e(amn)r](mn/x) 

m<D n 

gcd(m,tu)=l gcd(m,to)=l 

(3.6) 

= fi(m) fi(d) e(amdri)r)(mdn/a 

m<D d\w n 

gcd(m,ui)=l 



Let M = min(L-Dj , Q/2w). We start by separating the terms with m < M, 
divisible by q. Then e(amdn) equals e((8md/x)n). By Poisson summation, 



'^2e(amdn)r](mdn/x) = /( 



n 



where /(u) = e((5md/x)u)r]((md/x)u). Now 



/(n) = / e(—un)f(u)du = — - / e ( ( 6 ) u] rj(u)du = rj ( — -n — 5 

J md J \\ md) / md \md 



By assumption, md < Mw < Q/2 < x/2\5q\, and so \x/md\ > 2\5q\ > 25. Thus, 
by C3t, 



(3.7) = if 5M) + o,f E 1 



x ~, i- \ md In" 
—iVi-S) + 



2 °- Lr^i/2^ (n-r) 2 j 



mci a; (2ir) 



Since x H- 1/x 2 is convex on E+ , maX| r |<i/ 2 = Er^o (n-i/2)^ = 

7T 2 - 4. 

Write q = q + q~, where q + = gcd(q,w) and q~ = q/q + . We can assume q is 
square-free at this stage, since otherwise all terms with dm divisible by q vanish; 
then gcd{q + ,q~) = 1. Summing over all m < M, d\w with q\dm, we obtain, 
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using (|3.7p . that 



fj,(m) 5^ e(amdn)r](mdn/ x 

m<M d\w n 

gcd(m,ru)=l g|dm 

/x(m) x - /x(<i) 



y 

gcd(ro,w)=l g +|c( 
q~ \m 



+ o< 



\ 



i i \ |V|i 



4 7T' 



X, 



V 



m<M d\w 
q~\m q +\d J 



q w/q + 



] -rf(-5) ■ Yl 



n(m) 



m<M/q- 
gcd(m,wq~)=l 



( M(M + q-)q+a{w/q+) (1 __ _L_\ „ 
+ - V xq~ ' \8 2vr 2 J " 



»7 1 



Now consider the terms with md not divisible by q, or m > M. Here we apply 
2.ip with /(y) = rj(my/x) and am instead of a, and obtain 



(3.8) 



e(amn)r](mn/ x) 



< - 



77 1 1 • mm l — 



m sin i am 



where c = 2|r/|i/|r/|i + 1. 

Take first the case Q/2w > [D\; then M = \_D\. In this case, all terms with 
md divisible by q have been removed. The contribution of the other terms is 
bounded from above by 



1 1 /1 • 
-\rj \i ■ 2^ mm 

m<Dw 
q\m 



CX 



1 



m sin iram 



By Lemma 13.21 this is at most 

00 

E Yl min 

i = j q<m<min((j +l)q—l,Dw) 

\Dw/q]-l 



CX 



m ' I sin(7ram) 



<— log2eg + V* — log < — log 2eq + / -log^^dt. 



2q . lQcx 2q 
— log — — < — 

TT jq TT 



2 , 16cx 



Let us now examine the case Q/2w < [D\; then M = Q/2w. The sum of the 
terms in (j3.6j) with md not divisible by q is bounded from above by 



l \ 'I 



E 

m<M 



mm 



cx 



m sin fam 



+ E 



mm 



M<m<Dw 



CX 



m ' I sin(7ram) 
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As above, Lemma 13.21 gives us that the sum over m < M, q \ m is at most 



2q f M 2 16cx 



— log 2eq + I — log dt. 

TT J TT t 

Let us now treat the terms with M < m < Dw using the approximation a' jq' 
instead of a/q. Lemma [3TT1 gives us 

E(cx 1 
mm — , — - — - 
V m sm(TTam) 

M<m<Dw ' 
q\m 



<Y^ mm(- X ' 



j=0 jq'+M<m<m'm((j+l)q'+M,Dw) 
\(Dw-M)/q']-l 



\ m ' | sin(7ram) 



Ei 3cx 2q' 2ecx 



s [M + jq'\ +1 tt & [M+jq'\ +1, 

The first kind of terms within the sum makes a small contribution: since q' > 

(e/(l + e))Q and M = Q/2w, 

(3.9) 

E )/q ^ 3cx 3cx 1 f Dw 3cx 6wcx 3(1 + e)cx Dw 
^ ~M~+~Jq ! ~M + ~q' Jm ~ * = ~Q~ + ~^Q~ ^ ~W 

The second kind of terms add up to 

r(^-M)/ 3 H-i 2ecx 2q' Aewcx 2 f Dw , 2ecx J 

> — log : < — log — 1 — / log at. 

„ s L M + jVJ+l-vr S Q + 1 ttJ m t 

Recall that q' < Q' < (1 + e)Q + 1 < (1 + e)(Q + 1). Recall as well that Q/2 < 
w[D\; since Q is an integer, it follows that Q + l < 2wD. Since t \-+ tlog(Aecx/t) 
is increasing for t G (0, 4cx] and (Q + 1) < 2wD < 2x < Acx, 

. . 2o' Aewcx AwCl + e)D Aewcx Aw(l + e)D 4ecx 
(3.10) — log— — < log— — - = lo g^TT- 

TT Q + l TT 2wD TT 2D 

Summing all terms with m not divisible by q, we obtain at most 

2« 2 f Dw , 16cx , Wl + eYD, 2ecx 6wct 3(l + e)cs, Dto 

-Mog2eg+- / log-— cft+^ — log — — — I — — — h — -p — log — — 

7T TT J t TT D Q eQ M 

where the last three terms are present only if Q /2w < \_D\ . 
Now 

rDw , 16cx , 2 „ , 16ecx 
log at = —Dw log ■ 



2 r 

K Jo 



TT Dw 

We remember to multiply by (l/2)|r/|i (as in (]3.8[) ) and obtain 

(2 + 2e)Dw, 2ecx Dw , I6ecx\ 
log 2eg H log — — H log — — 

L/ 7T l«iJ / 



7T TT 



, /1N (3wcx 3(1 + e)x , DuA 

+ (2|ml + KI)-(— + ^^108—] 

with the slightly better bound 



1,11, 1 ^ 16ecx\ 

77! • - log 2eq + -Dw log — — 

1 TT TT Dw J 
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if Q/2w >[D\. □ 

We will also need a version of Lemma 13.31 with a factor of (log n) . 

Lemma 3.4. Let a = a/q+5/x, gcd(a, q) = 1, \5/x\ < 1/q 2 . Letrj be continuous, 
piecewise C 2 and compactly supported, with n" G L\. Let w G Z + be square-free. 
Assume that, for any p > po, po a constant, the function i]i :P (t) = log(pt)r]o(t) 
satisfies 

(3.11) \m, P \i < log(p)|r?|i, \v[Ji < log^k'li, WIJl < log(p)|»/'|i. 

Then, for D < xj max(ew, po, 4.33) and any e G (0, 1], 
(3.12) 

p(m) (log n)e(amn)r](mn/x) 



m<D 
gcd(m,«;)=l 



gcd(n,M))=l 



xp(q) 4>(w/q~ { 



w/q^ 



-?(-<) 



E 



+ 



x/i(g) 



(w/q + ) 



m<M/q- 
<rcd(m,wq~)=l 

( 



Mm) lo /^\ 



log-77(-<5) 



7T 



w/q 

log 2e<7 log x ) |t/'h 



E 



p(d) log <i 



7?(-<5) 



£ 

m<M/q~ 
gcd(m,iu<j — )=1 



p(m) 



2t(w) , ,D 2 w, 

——^+a(w) log 

e xg 



e^ 2 x 
D 



1 1 

8 ~ 2^2 



IV'Ii 



„ / / 3wcx 3(1 + e)cx , Dw\ / x \ . ,. \ 
+ ~W~ " "M J ( l0g m) 11/11 J 



3 + 2e 4ecx ex 
log — — log — 

7T D D 



where c = 2|j7| 1 /|t/| 1 + 1, Q = [x/\5q\\, M = mm([D\,Q/2w), q + = gcd(g, w) 
and q~ = q/q + ■ If Q/2w > D, the last two lines of 13.12\) are replaced by 



(3.13) 



W\i t — \ ( i 16cex . ex , 
LUlDw log — — log — + 1 . 
7T V L>w Dw 



Proof. Define Q, Q' , a', q' and M as in the proof of Lemma 13.31 The same proof 
works as for Lemma 13.31 we go over the differences. When applying Poisson 
summation or (|2.ip . use fjxxlmdip) = 0°S xt/md)rj(t) instead of n(t). Then use 
(|3.11|) to bound the l\ norms of r)i x / md and its derivatives; for example, 

\< x/md h<^ d \v"h. 

Then, for f(u) = e((Sm/x)u) (log u)r]((m/x)u), 

/(„) = r?1;c/m(i (- ( 5) + — 2-L_—0<(k 2 - 4) 



x (2ir)< 



^(,^ M , + 1 „ g (^) 5M ,) + ^,„ g (^) 



i»"ii i»"ii 



7T- 1 
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The part of the main term involving log(x/md) becomes 

m<M d\w 
gcd(m,io)=l q +\d 
q~\m 

= xn{q)(j){w/q+) _ fi(m) / x_ 

q w/q + ' m Vmg 

m<M/q~ 
gcd(m,wq~)=l 

-^slr}(-5). y aw y ^±io g d 

q ^ m ^— ' d 

m<M/q~ d\w/q+ 
gcd(m,wq~)=l 

It remains to see how the same factor affects the error terms. 
In the term in front of (|t?"|i/4 — \r]"\i/ir 2 ), we find the sum 

, / x \ fx f M/q ~ , x/(q-d) \ x M 2 e^x 

X>M^) ^ (^p+jf tlog-^pcftj = s + — log — , 

where we use the fact that i h-> ilog(x/t) is increasing for t < x/e. By the same 
fact (and by M < D), (M 2 /q~) \og{e 1 ' 2 x / dM) < {D 2 /q~) \og{e 1 ' 2 x / dD) . 
The basic estimate for the rest of the proof (replacing (I3.8jl ) is 



y^e(amn) (log n)q (j^) =^2e(amn)r] ltX/m (™J 



1 1 / 1 ( cx 1 

< ~\ r ] ll x/m ' mm > 1 ? s 

2 ' ' \m | sin(-7ram) 

log x/m . ( cx 1 

< 1?7 |i • mm 1 



m sin iram 



The term (g/7r) log 2eg becomes (q/n) log x log 2eq; the terms in (|3.9p get mul- 
tiplied by logx/M. Recall that (Q + 1) < 2wD. The analogue of (l3~10jl is 

2g' 4euicx , 2wx 4(1 + e)wD , 4ecx 2x 

~V g qTT g qTT n log ^F log 2D' 

which holds because D < x/4.33 and c > 1 imply that we are in the range where 
the function t h-> f log(4eu;cx/i) log(2u;x/t) is increasing. The other term of this 
type is the one coming from the integral 

2 f Dw , 16cx , a; , 2 „ / 16cex , ex 
- / log — — log -dt = -Dw log — — log — — + 1 

TT J t t 7T \ LJW Vw 

The two terms, summed together, give less than 

6 + 4e 4ecx ex 
Dwlog— — log— , 

7T D D 

which is then multiplied by (l/2)|?/|i. □ 
We will apply the following only for q relatively large. 
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Lemma 3.5. Let a = a/q + 0<(l/q 2 ), gcd(a,q) = 1, q > 2e 4 . Let rj be con- 
tinuous, piecewise C 2 and compactly supported, with rj' € L\. Let w G Z + be 
square-free. Then, for any r > 1, U, V > 1 with UV < x/w and e < 1, 



(3.14) 



A(w) /u(n) e(avun)rj(vun/x) 

v<V u<U n 

gcd(«,w)=l gcd(M,to)=l gcd(n,u>)=l 

< %mi*<^ iog(7 S -) + ^ - 4) . i|qL^/' H 

q w/q + (2tt) z x 



1.03883 



w/q+ ~ ov ' s ' ' v " ^ (2vr) 2 

t/y 2 (iogy)(?7y) 2 u 2 + u 



+ 



2q- 



(logq )q~ 



1, /, /^2<? , , 2 / C/F , 16ecx 
+ « *7 i • — log 2eg log g + -C/F log log — — - 1 

Z \ 7T tt \ e U V 

1 , Aq 
+ « 1*7 |i - — log{UV + q) log 16cx, 

Z 7T 



i/E7V < Q/2tw, where c = 2\q\x/\rf\ 1 + 1, Q = |x/NJ, ^ 
g - = g / 9 +. 7/C/y > Q/2w, 
(3.15) 



A(i>) / i ( n ) e(avun)r](vun/x) 

v<V u<U n 

gcd(«,«))=l gcd(u,u>)=l gcd(n,uj)=l 

x|r?((5)| 4>(w/q + ) 



icd(q,w) and 



< 



w/q^ 



\og(Vq-) 



1. „ 2 / C/y™ 16ecx 1 

+ o *7 i • -*7F«; log log — — - 1 

Z 7T \ e (/Kio 

14 / 2ecx \ 

H — l*7'|i • — (1 + e)UVw I log UVw ■ log + log 2ecx • max(logu>, 1) I 
2 7r \ c/V / 

+ nW\i ■ \—^og(UV + g) log 16cx + — log2e g logg 

2 \ TT TT 

,2 ^ Irfh <r(w) f QV logV ( Q \ 2 U 2 + U n 

+ ^ |*/|l • + log j log 



Proof. We proceed essentially as in Lemma 13^3] let Q, q' , Q', M, t/ + and g be 
as there. The sum can be written as 



A(f) ^2 ^{u)' S ^ l ^(d)' S ^ l e{avudn)r](vudn / x) . 



v<V 
gcd(i>,«j)=l 



u<[7 
gcd(u,«j)=l 
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Let M = mm([UV\ , Q/2w). The terms with uv < M, uvd divisible by q 
contribute 
(3.16) 

/ \ 



E 

a<M 



E A (v)v(a/v) 



v\a 



gcd(a,w)=l \ a /U<v<V 
q-\a 



d\w 
I q+\d 



x v( s ) , n ( ad , 2 



x (2?r) 



The contribution of xrj{5) / ad can be rewritten as 
(3.17) 

/ 



xrj{S) 



+ xrj(5) 



4>{w/q + ) 
w/q + 

<l>(w/q+)n(q+) 1 



E 



A(v) 



E 



H{u) 



\ 



v 1 — ' u 

v<V u<mm(U/g- ,M/Vq~) 

gcd(v,wq~)=l gcd(u,wq~)=l 



w/q^ 



p a <V 

p a \q~ 



■ - E a* (logp) 



E 



u<min(U/(q- / p a ) ,M / (V q~ / p a )) 
gcd{u,wq~ /p a )=l 



U 



We are happy to bound the innermost sum by 1 (using (|2,5p ). and obtain a total 
of at most 



(3.18) 



x\rj(5)\ (j){w/q + ) 
q w/q + 



(log V + log q 



by (I2TTID : recall also that rj(6) < min(|??|i, \t]"\ 1 /(2tt5) 2 ) (see QEty). 
The second term in f|3. 16j) contributes 

(^-&: EE A(tOj>d 



(2vr) 2 x 



(vr 2 - 4) 



u<Uv<V d\w 

gcd(uv,w)=l q +\d 
uv<M, q~\uv 
u sq-frcc 

r/'\i q + a(w/q + ) 



(2vr) 2 x 



E E A ( v ) uv - 

u<Uv<V 
gcd(uv,w)=l 
uv<M, q~\uv 
u sq-frcc 



By the estimates (|2.10p and (|2.12p , the double sum here is at most 
(3.19) 

Y A ( ,-)v mhl{U > M ^)( min (^ M/v) + q~) + y U(U + q-/p) 

^ 2 v~ p% 2 <i lv 



v<V 
gcd(n,g-)=l 



p\q~ 



< 1.03883^— + -^log(T/)M 2 + - ±£ (logg~)g- 



MV 1 
1~ + 2q~- 



(A stronger bound, without the first term, holds for q non-square- free.) 
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From this point onwards, we use the easy bound 
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A(u)/i(a/t 

v\a 
a/U<v<V 



< log a. 



What we must bound now is 

(logm) ^2 e{amn)rj{mn / x) . 



m<UV 
gcd(m,w)=l 
q j m or m > M 



gcd(n,«j)=l 



Now use (|3.8p and Lemmas 13. 1H3.2I The terms with m < M, q { m contribute 
at most 



(3.20) 



*l 'I 



29 1 , r ^ ^ 2 9 , 16cx 
— log2eglogg + > (log(j + l)g) log — 

7T ^— ^ 7T J<7 



For any i G [a, 6] , < a < b < y, 
(3.21) 



y y 

log b log log t log - 

a t 



log - log — + log t log — 

t a a 



y 



< log max I — , t) • log — < log y • log - . 

Hence (taking out the term with j = \M/q] — 1) we see that the sum in (|3.20p 

is at most 

(3.22) 



M n , 16c* 2« 16cx 2g rM ^I\ j + 1 _ 
log t log est H log(M + g) log 1 > log log locx 

t IT q 7T frf j 



2M 



< — I log — lo. 

7r V e 



M , 16ecx 



M 



4q 



1 + — log(M + q) log 16cx. 



7T 



If [f^^J 5; Q/2w, then we have considered all terms, and get a total of at most 

%*>l^£W<r) + (. 2 - 4) ■ 

q w/q + (2tt) z x 



1.03883 



^ 1 v " (2vr) 2 

in/ 2 (io g y)(?7y) 2 u 2 + u 



2 

7T 



+ 



2q 



Qogq )< 



log 2eg log q H £/V log log 



7T 



C/V , 16ecx 



uv 



1 4</ 

H — l^'li " — \og(UV + q) log 16cx, 

2 7T 

since, in that case, M = [UV\ . 

Assume from now on that [UV\ > Q/2w. Then M = Q/2w. The terms 
with m < M, q\ m contribute (|3.20p as before, whereas the terms with m > M 
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contribute at most 



(3.23) 2 1 1,1 ^ V L^J + + 1 L M J + h' + 1 

• logmin([Mj + (j + l)g',Dw), 

where D = UV. As in the proof of Lemma [373| q' > (e/(l + e))Q. 
The first kind of term within the sum in (|3.23p adds up to at most 

3c* /6wc* 3(l + e)c* Dw 

3.24 logDw- 2^ <l og Dw [——+ ' log — 



much as in f|3.9[) . We wish to approximate 

2g^- 1 / , 16c* 
vr ^ jg 

(3.25) 

r(^-M)/ 9 'i-i 2ecx 

+ V E (logmin(LMj + (j + lV,Dw))log [MJ - jg/ - i 
by an integral. By (13,2ip . the first sum is (as in (|3.22p ) at most 



2 

7T 



f M 16c* 4g 

/ log t log dt H log(M + g) log 16c*. 

Jo t 



The second sum, again by (|3.2ip . is at most 



Dw , , 2ec* , 2g' , „ , 2ec* 
log t log dt H log Dw log ■ 



7T J M t IT ° ° [M\ + 1 

2 „< rc^-wi-. [MJ + (J + ltf + 1 

+ V log i i ry— : log2ec*. 



j=0 

Clearly 

f M , , 16c* , /- Du ' , , 2ec* , f Dw , i 16c* 

/ log r log at + / log t log at = / log t log 

Jo * /m * Jo t 



dt 



. , 16ec* 
< Dw log log 1 

e Dw 



We conclude that (13.25P is at most 



2 „ /, flw, 16ec* \ 4g , 

—Dw log log — 1 J H log(M + gj log 16c* 

(3 26) vr V e / vr 

1 ' J 2g' 2ec* 2g' Dw 

+ - log Dw log + - log log 2ec*. 

Note that (much as in the proof of Lemma [373]) q' < (1 + e)(Q + 1), Q < 2[Djw 
(implying Q + 1 < 2[D\w < 2Dw) and M = Q/2w (implying [M\ + 1 = 
[Q/2w\ + 1 > (Q + l)/2w and M < D = UV). Proceeding as before (I5TTJ1) . we 
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obtain that ()3.26j) is at most 



2 / 16ecx \ Aq 

—Dw log log — 1 M logiUV + q) log locx 

ir \ e Dw J 7T 

ADw . ( 2ecx 

H (1 + e) log Dw log — — — h log 2ecx max(log w,l) 

7T \ D 



□ 



4. Type II 



We must now consider the sum 



(4.1) 



E 

m>U 



E^ d ) 



^] A(n)e(amn)r](mn/x). 

gcd(m,ui)=l \d|r 

Here the main improvements over classical treatments are as follows: 

(1) obtaining cancellation in the term /x>[/ * 1, leading to a gain of a factor 
of log; 

(2) using a large sieve for primes, getting rid of a further log; 

(3) exploiting, via a non-conventional application of the principle of the large 
sieve (Lemma 14. 3p . the fact that a is in the tail of an interval (when that 
is the case). 

Some of the techniques developed for ([T} should be applicable to other instances 
of Vaughan's identity in the literature. 

It is technically helpful to express rj as the (multiplicative) convolution of two 
functions of compact support: 

(4.2) n(x) = / m (t) m (x/t)-. 

Jo 1 

For the smoothing function rj(t) = rjo(t) = 4max(log2 — | log 2t| , 0) , (|4.2p holds 
with 



(4.3) 



2 if t€ (1/2,1] 
otherwise. 



We will work with rji(t), T]2if) as in (14. 3p for convenience, yet what follows should 
carry over to other (non-negative) choices of rji and rj2 ■ 

By P~2j) . the sum (|4TT|) equals 
(4.4) 

/ \ 



E 



m>U 



n>V 



E^ d ) 

... . ■ d>U 
gcd(m,w)=l \d\m / gcd(n,u))=l 

/ 



x/U 



A(n)e(amn)r)i(t)rj2 
\ 



mn/x\ dt 

t / T 



V 



A(n)e(amn )— — 
W 



E E^) E 

ma, x (^,U)<m<^ \^> U J max(V,-f)<n<VK 

gcd(m,«))=l < ^ m ^ gcd(n,w)=l 

by the substitution t = (m/x)W. (We can assume V < W < x/U because 
otherwise one of the sums in (|4.5p is empty.) 
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We separate n prime and n non-prime. By Cauchy-Schwarz, the expres- 
sion within the integral in (|4.4p is then at most y/S\JU, W) ■ S2QJ, V, W) + 
y/Si(U,W) ■ S 3 (W), where 

/ \ 2 



Si(U,W) 



E 



(4.5) 



<(iwV)<™<w \ d > u 

. . \d\rn 
gcd(m,ui)=l 



/ 



S 2 (U,V,W) 



E 



gcd(m,u))=l 



E 



(log p)e (amp) 



max(V,-^)<p<W 
gcd(p,?o)=l 



and 



(4.6) 



E 

gcd(m,tu)=l 

E 



E A ' n ; 



n<VK 
n non-prime 



1.42620W 1/2 ) < 1.0171a; + 2.0341VF 



gcd(m,iy)=l 

(by [RS621 Thm. 13]). We will assume V < w\ thus the condition gcd(p, w) = 1 
will be fulfilled automatically and can be removed. 

The contribution of S^(W) will be negligible. We must bound Si(U,W) and 
S 2 (U,V,W) from above. 

4.1. The sum Si', cancellation. We shall bound 



Si(U,W) 



E 

max(U,x/2W)<m<x/W 
gcd(m,w)=l 



E ^ 



d>U 
d\m 



J 



There will be what is perhaps a surprising amount of cancellation: the expres- 
sion within the sum will be bounded by a constant on average. 

4.1.1. Reduction to a sum with fi. We can write 
(4.7) 



E 



\ 



E^) 



m&x(U,x/2W)<m<x/W 1 d>U 
\d\m 



E E Kdi>U)^d 2 >U) 



,d2\m 



gcd(m,w)=l 



gcd(m,io)=l 



E E E M(nOMM E !. 



ri <x/WUr 2 <x/WU I 

gcd(n,r 2 )=l gcd(/,nr 2 )=l 
gcd(nf-2,m)=l nl,r 2 l>U 
gcd(£,w)=l 



X ^ ^ X 

2W < - m ^ w 

r\T2l\vn 
gcd(m,u))=l 
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where we write d\ = r±l, di = r 2 Z, I = gcd(di, d^). (The inequality r\ < x/WU 
comes from ri^lm, m < x/W, r%l > V r ; Ti < x/WU is proven in the same way.) 
Now (077D equals 



(4.8) 



Y Y Y / x ( r i)/ i ( r 2) 



Y 



S <W77 '-1<W777 ^2< 



WUs A WUs 



gcd(s,ui)=l gcd(ri,r 2 )=l 
gcd(rir 2 ,iu)=l 



max! 



U 



x/W ) <1< x/W 



gcd(Z,r 1 r 2 )=l, (^(/)) 2 =l 
gcd(£,u>)=l 



where we have set s = m/{r\riV)- 
Lemma 4.1. Let z,y > 0. Then 



(4.9) 



ri<y r2<y 
gcd(ri,r 2 )=l 
gcd(rir 2 ,ui)=l 



5 z 

V min(r^ ,r2) ' 2r^ r 2 / — r l r 2 

gcd(«,r 1 r 2 )=l,( At (0) 2 =l 
gcd(£,iu)=l 



equals 



(4.10) 



6z w ^ ^ MOlM 7 ^) 



-, / / / \ / \ I lUcLX I 

Tro-(w) £< v ££ y °"(n)CT(r 2 ) V V 2 y y 



gcd(n,r 2 )=l 
gcd(rir 2 ,«j)=l 



+ 0< 5.08 C 



If w = 2, the error term in ^4-10 ) can be replaced by 



(4.11) 



o s u. J 7<['|) , »vi.(i + -^) (i-^; 2 



Proof. By Mobius inversion, (|4.9p equals 
(4.12) 

ri<» r 2 <?/ K 
gcd(ri,r 2 )=l 
gcd(rir 2 ,u!)=l 



£>min 



r l r 2 

z/y 



min(r^,r 2 ) ' 2r^r 2 

;cd(l,w)=l 



Y Kdi)Kd2) 

d 1 \r 1 ,d 2 \r2 
did,2\l 

Y^ d3 "> Y ^ m )- 

d3\w m 2 \l 

d$\l gcd(m,rir 2 «j)=l 



We can change the order of summation of rj and di by defining = ri/di, and 
we can also use the obvious fact that the number of integers in an interval (a, b] 
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divisible by d is (b - a)/d + 0<(1). Thus (j4TT2l) equals 
(4.13) 

E /x(di)/x(d 2 ) E n(d 1 si)n(d 2 s 2 ) 

di,d 2 <y s 1 <y/di 
gcd(di,d2)=l s 2 <y/d 2 
gcd(d 1 d 2 ,w)=l gcd(disi,<i2S2)=l 

gcd(sis 2 ,«))=l 



E /*wo 

d,3\w m < 



£ 



n(m) 



d'^s 1 d'^s 2 d' i 



d\d 2 d 3 m 2 s\d\s 2 d 2 



1 — max 



1 s\d\ s 2 d 2 

2 y y 



plus 



(4.14) 



gcd(m,c(isid2S 2 ui)=l 



0< 



\ 



E 

<2i,d 2 <?/ 

;cd(did2,iu)=l s 2 <y/d 2 

gcd(sis 2 ,w)=l 



£ £ £ 



1 



si<j//di d 3 |«j m < 



I Pl a 2 s 2 a 3 



m sq-free y 

If we complete the innermost sum in (|4.13p by removing the condition m < 
z/ (df sd^), we obtain (reintroducing the variables rj = djSi) 

n7?< y rir2 v V2 y y 

gcd(n,r 2 )=l 
(4.15) gcd(rir 2 ,ui)=l 



EE E 

^2^2 



fi(di)n(d 2 )n(m)fi(d 3 ) 
d\d 2 d^m 2 



]cd(rn,rir 2 w)=l 



times 2. Now (|4.15|) equals 
fJ,(ri)(J,(r 2 )z 



E 



n,r 2 <y 
gcd(n,T , 2)=l 
gcd(rir 2 ,ui)=l 



nr 2 



6z to 



1 — max 



1 T\ r 2 

2 y y 



n K)nH 



p|rir2W 



p\rir 2 



Li(ri)n(r 2 ) / An r 2 

> / — 7 — \ — 7 — 7 I ITlctX I — , — ; — 

^ ri j? <y °"(n)o-(r 2 ) V V 2 y y 



E 



gcd(ri,r 2 )=l 
gcd(rir2,ui)=l 



i.e., the main term in (|4.10p . It remains to estimate the terms used to complete 
the sum; their total is, by definition, given exactly by (|4.13p with the inequality 

m < \J ' z I \d\sd\s<id%) changed to m > y/z/(d'(sd'is 2 d 3 ). This is a total of size at 
most 



(4.16) 



E 

dx,d 2 <y 



E E E 



si<y/di d 3 \wmy / z ~ 
gcd(d 1 d 2 ,w)=l s 2 <y/d 2 V ■<-rv'i~rh 

gcd(s 1 s 2 ,w)=l 



d\d 2 d 3 m 2 s\dis 2 d 2 



m sq-free 
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Adding this to (|4.14p . we obtain, as our total error term, 
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(4.17) 



where 



E EE/ 

d\,d 2 <y si<y/d\ ds\w 

gcd{d 1 d 2 ,w)=l s 2 <y/d 2 

gcd(sis 2 ,«i)=l 



d\ sidls 2 d 3 J ' 



/<*> = E 1+5 E i 



m sq-free 



2 £ — ' 

m>x 

m sq-free 



It is easy to see that f(x)/x has a local maximum exactly when x is a square-free 
(positive) integer. We can hence check that 



/(•<•; < ~ ( - + - } 



02) 
C4 



1.25 \ \x = 1.26981... x 



for all x > by checking all integers smaller than a constant and using {m : 
m sq-free} C {m : 4 { m} and 1.5 • (3/4) < 1.26981 to bound / from below for x 
larger than a constant. Therefore, (j4.17j) is at most 



^ E EE 

d\,d 2 <y si<y/di ds\w 

gcd(did 2 ,w)=l s 2 <y/d 2 

gcd(sis 2 ,w)=l 



d\sid\s2d% 



( 



p\w 

We can bound the double sum simply by 



\ 



S S d^ 

d<y s<y/d 
\gcd(d,to)=l gcd(s,io)=l / 



E E^E^v^)n(i-^). 

d<y s<y/d V d<y v 7 p\w V L 7 



<2<y s<y/d 
gcd(d,w)=l 

Alternatively, if tu = 2, we bound 



1 1 

E Ts = E vs sl + 2 

gcd(s,w)=l s odd 



1 1 



ds = yy/d 



and thus 



S Yl ( ] 

d<y s<y/d d<y 

gcd(d,w)=l gc d(s,w)=l gcd(d,2)=l 



23/2 ) C ( 2 ' ' 



□ 
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Applying Lemma [4.11 with y = S/s and z = x/Ws, where S = x/WU, we 
obtain that (I4.8P equals 
(4.18) 

6x w \ - 1 \ - \ - ^(ri)/x(r 2 ) / / 1 _n ^2_ 

V ; s<S' n<S/sr 2 <S/s V -L/ V -V \ \ / / , 

gcd(s,TO)=l g cd(ri,r 2 )=l 
gcd(rir2,tf)=l 



\ p\w v 



1 - — X 

p3/2 



with 5.04 replaced by 1.27 if w = 2. The main term in (|4.18p can be written as 

( 4 - 19 ) E 1 T E E 4^4^- 

^ s ^ s A/2 ^ ^ a(n)a(ra) 

gcd( S ,w)=l gcd(ri,r 2 )=l 
gcd(rir2,i«)=l 

From now on, we will focus on the cases w = 1 and w = 2 for simplicity. (Higher 
values of u; do not seem to be really profitable in the last analysis.) 



4.1.2. Explicit bounds for a sum with fj,. We must estimate the expression within 
parentheses in (|4.19p . It is not too hard to show that it tends to 0; the first 
part of the proof of Lemma S2] will reduce this to the fact that E^ n //(n)/n = 0. 
Obtaining good bounds is a more delicate matter. For our purposes, we will need 
the expression to converge to at least as fast as l/(log) 2 , with a good constant 
in front. For this task, the bound (12.6P on ^2 n<x [i(n)jn is enough. 

Lemma 4.2. Let 

«.(«)=- EE ! Wr 1 v 

o-{ri)a(r 2 ) 

ri<x r2<x 

gcd(n,r2)=l 
gcd(rir2 ,w)=l 

where w = 1 or w = 2. Then 



\9i(x)\ < < 



\92(x)\ < < 



'l/x if33<x<10 6 , 

i(111.536 + 55.768 log x) ifW 6 < x < 10 10 ; 

2^0044325 + 0T079 2 / x >10 10 , 

v (logx) 2 y/x J — ' 

2.1/ x i/l^x^lO 6 , 

±(1634.34 + 817.168 log x) if 10 6 <x< 10 10 ; 

00381^ + 02046, i/x>10 10 . 

(log 2+ yhc J — 



The proof involves what may be called a version of Rankin's trick, using Dirich- 
let series and the behavior of £(s) near s = 1. The statements for x < 10 6 are 
proven by direct computation. 



MINOR ARCS FOR GOLDBACH'S PROBLEM 



29 



Proof. Clearly 



(4.20) 



EE E 

ri<xr2<x \d\ gcd(n,r2) 
gcd(rir2,ui)=l 

= E m EE 

d<x ri<x r2<x 

gcd(d,w)=l d\ gcd(ri,r-2) 

gcd(rir2,«))=l 

M (d) 



/f(ri)/x(r 2 ) 
cr(ri)o-(r 2 ) 

cr(ri)o-(r 2 ) 



\ - \ - //(m)/x(M 2 ) 

(<r(d)) 2 ^, ^, a( Ul )a(u 2 ) 

d<x y y " Ul <x/d u 2 <x/d V ' V ' 



E 

<i<a; 
gcd(d,w)=l 



gcd(ui,dw)=l gcd(u2,dw)=l 

I v- 



E 

d<x 
gcd(d,w)=l 



E 



H(r) 



a r 



r<x/d 
\gcd(r,d«i)=l / 



Moreover, 



E 



r<x/d 
gcd(r,dui)=l 



o-(r) 



y- A t ( r ) v^- tt / P 

r<:r/(2 ci'|r p|<2' V ^ 

gcd(r,dw)=l 



d'<x/d \p\d 
fi(d') 2 =l 
gcd(d' ,dw)=l 

- E 1 



- e. (n^) e 

r<x/d 
gcd(r,dw)=l 
d'\r 



E 



d'<x/d v ' r<x/dd' 

M (d')2=i gcd(r,dd'w)=l 
gcd(d',dui)=l 



and 



Hence 
(4.21) 



E 



r<x/dd' 
gcd(r,dd'w)=l 



E ^ E 

d"<x/dd> r<x/dd'd" 
d"\(dd'w)°° 



E $f 

gcd(d,ui)=l 



E ^ E 

d'<x/d V y d"<x/dd' 

fj,(d') 2 =l d"\(dd'w)°° 
\gcd(d' ,dw)=l / 



where /(i) = £ r < t MMA • 

We intend to bound the function fit) by a linear combination of terms of the 
form , <5 G [0, 1/2). Thus it makes sense now to estimate F w (s\, s 2 , x), defined 
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to be the quantity 



E 



(MM 

((7(d)) 5 



gcd(d,w)=l 



E If, E ^w«' 



-si 



\gcd(c( / 1 ,d«i)=l 
/ 



di'|(dditu)°° 



E 



\gcd(c(2i^)=l 

for si, S2 £ [1/2, 1]. This is equal to 



2 v 27 ■ 



E i • 



'/ J//\1-S2 



E 



« '""'Id+p-'fd-p-'in, 

gcd(d,«i)=l 



i 



E 
E 



'A 2 



|w (l-p- s l)(l-p- s 2) 
\ 



(1-p- 



-*2 



n 



i 



w i+l ^ l j(i+p'- l )(i-p'- 1 ; 



'A2 



n 



i 



(^•2+1 11 (1 +J /-1) (1-^-2 

\gcd(d',cto)=l 



/ 
\ 

I 



which in turn can easily be seen to equal 

nfi+ p -' lp " a in 



(4.22) 



(1 + J> _1 )(1 -P" S2 + P' 1 ) J 11 (1 -P~ S1 )(1 -P~ S2 ) 



pfu> 



n (i+ 



p 1 p si 



(l +p -i)(l_p-.i; 



n i+ 



J) Ij, s 2 



(1 +P~ 1 )(l ~P~ S 



S2 



Now, for any < x < y < x 1 / 2 < 1, 

(l+x-y)(l-xy)(l-xy 2 )-(l+x)(l-y)(l-x 3 ) = (x-y)(y 2 -x)(xy-x-l)x < 0, 

and so 
(4.23) 



1 + 



xy 



(l + x-y)(l -xy)(l -xy 2 ) 



< 



(1 - x 6 ) 



(1 + x)(l - y) (1 + x)(l - y)(l - xy)(l - xy 2 ) (1 - xy)(l - xy 2 ) ' 
For any x < yi, y 2 < 1 with y 2 < x, y| < x, 

(424) 1 + ™* < (l-x 3 ) 2 (l-x 4 ) 

1 " (1 - yi + x)(l - y 2 + x) - (1 - y iy2 )(l - y iy 2 )(l - y\y 2 ) ' 

This can be checked as follows: multiplying by the denominators and changing 
variables to x, s = y± + y 2 and r = y±y 2 , we obtain an inequality where the left 
side, quadratic on s with positive leading coefficient, must be less than or equal 
to the right side, which is linear on s. The left side minus the right side can be 
maximal for given x, r only when s is maximal or minimal. This happens when 
Vi = V2 or when either yi = ^Jx or yi = x for at least one of i = 1, 2. In each of 
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these cases, we have reduced (|4.24p to an inequality in two variables that can be 
proven automatical!)! 2 ] by a quantifier-elimination program; the author has used 
QEPCAD [HBllj to do this. 
Hence F w (si,S2,x) is at most 

(l-p-^(l-p-4) 1 

J- J- (1 _p-81-32)(l -p-2 S1 -S2\n -pSl-2s 2 \ 11 H _ p -Sl\n -pS2) 

p\w p\w 

( 4 - 25 ^ n (i + p -si-i)(i + p -2si-i) n (i + p -«-i)(i + p-2 fl 2-i) 

r C(gi + l)C(^ + l)C(2si + l)C(2g 2 + l) 

C(3) 4 C(4)(C(si + S2)C{2si + S2 ) C(si + 2S2)) -1 > 

where 



Jl if to = 1, 

I (l-2- 3 l+"2)-l(l-2- 2s l-"2)-l(l-2- 3 l)(l-2- s 2)(l-2-») 4 (l-2- 4 ) 11 W ~ Z - 

For 1 < t < x, (H2D and imply 



if x < W 1 

— log log ID 10 

2 + oo| (xy ogx _ logWll) if a . > iQiO j 



(4.26) /(t) < 



where we are using the fact that log a; is convex-down. Note that, again by 
convexity, 

log log x — log log 10 10 ... 1 , , A 

log, -log IP < °° g " I — '»"' = togloio = 00434294 ■ ■ ■ 

Obviously, y/2/t in (|4.26[) can be replaced by (2/t) 1 / 2_e for any e > 0. 
By (IQTj) and gj26]), 

|^(x)|<(-J i^(l/2 + e,l/2 + e,x) 

for x < 10 10 . We set e = 1/logx and obtain from (|4.25p that 

+ , 1/2 + , rXT C(1 + 26)C(3/2) 4 C(2) 2 

(427) F(1/2 + e> 1/2 + e ' x) - — mw) — - 

< 55.768 • C w ^ +e ■ (l + ^ 
where we use the easy bound Q{s) < 1 + l/(s — 1) obtained by 

/oo 

(For sharper bounds, see |BR02j .) Now 

(1 - 2- 3 / 2 -<)' 2 (l + 2- 3 ' 2 ) 2 (l + 2- 2 ) 2 (l - 2- 1 - 2 ') 



2 In practice, the case yi = y/x leads to a polynomial of high degree, and quantifier elimination 
increases sharply in complexity as the degree increases; a stronger inequality of lower degree 
(with (1 - 3a: 3 ) instead of (1 - x 3 ) 2 (l ~ x 4 )) was given to QEPCAD to prove in this case. 
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whereas C, 1 , i, = 1. (We are assuming x > 10 6 , and so e < l/(logl0 6 ).) 
Hence 

f i (111.536 + 55.768 log x) if w = 1, 
[5(1634.34 + 817. 168 log a;) ifw = 2. 

for 10 6 < x < 10 10 . 

For general x, we must use the second bound in (|4.26p . Define c = l/(log 10 10 ). 
We see that, if x > 10 10 , 

0.03 2 

\g W (x)\ < -TT- ^1(1 - C, 1 - C) • C w ,1-c.l-C 

(log a;)^ 



+ 2 ^003 

7 x log X 



+ -(111.536 + 55.768 log x) ■ C„ 



For w = 1, this gives 



Un(.)| < T^f + ^ + i(H1.536 + 55.768 log,) 

(iOgXj^ y/XlOgX X 

„ 0.0044325 0.1079 
~~ (log x) 2 



for w = 2, we obtain 



, 0.038128 25.607 \, , 
52 x < r^ + ^ + - 1634.34 + 817.168 log z) 

(logx)" 1 vxlogx x 

„ 0.038128 0.2046 
~~ (logx) 2 sfx 



□ 



4.1.3. Estimating the triple sum. We will now be able to bound the triple sum 
in (j4"49j) . viz., 



1 f 1 

(4.28) V - / g w {uS/s)du 

s<S S Jl / 2 
gcd(s,«j)=l 

where g w is as in Lemma 14.21 

As we will soon see, Lemma l4.2l that (14.28H is bounded by a constant (essentially 

because the integral J^ 2 l/i(logi) 2 converges). We must give as good a constant 
as we can, since it will affect the largest term in the final result. 

Clearly g w (R) = g w ([R\). The contribution of each g w (m), 1 < m < S, to 
(|4.28p is exactly g w (m) times 
(4.29) 

Y r(m+l)s/S 1 r(m+l)s/S 



1 ,1 ! Am+X)s/S 

E 7 / J u+ E - / ,„ E 



1 +1 <S — m 2171^-111+1 2(m+l) <S — 2m 

gcd(s,to)=l gcd(s,ui)=l gcd(s,«i)=l 

1 m\ \ - 1 v - /m + 1 1 



s -'ms/S s " s a s '""V s «Jl/J 



E i-SH E i + E 



^8 SJ s ^ s S s ^ s \ S 2s 

- J tt< s< — T L < s< ~ S nr o? _^ <s<#- 

m+1 — m 2m — m+1 J(m+1) — 2m 

gcd(s,w)=l gcd(s,uj)=l gcd(s,iu)=l 
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Write f(t) = 1/S for S/2m < t < S/(m + 1), f(t) = for t > S/m or t < 
5/2(m+l), f(t) = l/t-m/S for S/(m+l) < t < S/m and f{t) = (m+l)/S-l/2t 
for S/2(m + 1) < t < S/2m; then (021 equals En: g cd(n,«,)=i /( n )- % Euler " 
Maclaurin (second order), 
(4.30) 

£./»= f{x)- l -B 2 {{x})f{x)dx= r f(x) + 0< (hf"( x )\\dx 



-oo 

fOO 



poo 1 
/ f(x)dx + - ■ 0< 

J-oo 6 



1 2m 



1 / 1\ 1 n [ 2m 

_ log 1+ + -.o< - 




Similarly, 



/OO 1 
/(2x + l)--Sa({x}) 
-oo ^ 



d 2 /(2x + 1) 



n odd 



' — oo 

2 



/(x)dx - 2 / -B 2 

oo J—oo 
;>oo ^ /-oo 



dx 2 

X — 1 



f{x)dx + 



6 



2 VI 2 

0< (|/"(x)|)dx 



/ (x)dx 



log 1 



my 3 _ \V S / 



We use these expressions for m < Co, where Co > 33 is a constant to be 
computed later; they will give us the main term. For m > Co, we use the bounds 
on |<?(m)| that Lemma 14.21 gives us. 

(Starting now and for the rest of the paper, we will focus on the cases w = 1, 
w = 2 when giving explicit computational estimates. All of our procedures would 
allow higher values of w as well, but, as will become clear much later, the gains 
from higher values of w are offset by losses and complications elsewhere.) 

Let us estimate (I4.28p . Let 



c w ,o 



1/6 iiw = l, 
1/3 if w = 2, 



1 if w = l, 
2.5 ifu> = 2, 



C w ,2 



55.768... 
817.168. . . 



if w = 1, 
if w = 2, 



111.536. 
1634.34. 



if w 
if w 



2. 



J 0.0044325 .. . if io = l, I 0.1079... if w = 1, 

CwA ~ 1 0.038128 ... if w = 2, Cu, ' B ~ | 0.2046 ... if w = 2. 
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Then (j4T28l) equals 



V- , n f<f>W, (, l \ r, ( 5m 2 + 2m + l\\ 



m<Co 



S/10 6 <s<S/C 1/2 



0< ( -^L ) du 



5/10 10 <s<5/10 6 



uS/s / 



s<5/10 10 

which is 



^ / 1 \ ^ / 5m 2 + 2m+l ^ 
2^9w{m)-—- log(l + - + 2^ I^MI-°< K.o e2 J 



m<Co m<Co 

^ / log 2 log 2, , , fi .v 2-\/2 

+ 0< [ ^l-qj" + "TTje" ( C ^,3 + C W , 2 (1 + log 10 )) + 1QlQ/2 C W , 5 



+ o<[ £ c "' 4/2 



„ s(log5/2s) 2 

s<5/10 10 V 6 ' ; 

for 5 > (Co + 1). Note that E,<s/ioi" s(io g s/2^ = Jo 2/1 ° 10 



Now 

c ui,4 



2/1010 1 dt _ c wA /2 { 0.00009923. . . if w = 1 

2 7 t(logt) 2 ~ log(10 10 /2) ~ j 0.000853636 ... if w = 2. 

and 

log 2 , . fi ., 2-\/2 f 0.0006506. . . if w = 1 

W (C - 3 + C - 2(1 + bg 10 )} + "Tb^ 5 = \o.009525 ... if w = 2. 

For C = 10000, 

1 / n , ( l^ f 0.362482... if «; = 1, 

— — - > 3toM • log n — = < 

«J 2 4^ V "V 1 0.360576... if w = 2, 

m<Co ' \ 

. , 2 , [6204066.5... if w = 1, 

Ci« o / 9wi ("i 5m +2m + l)<< 

' 4^ " I 15911340.1... if to = 2, 

m<Co K 

and 

'0.00006931... if 10 = 1, 



c^i • (Iog2)/C 



0.00017328... ifw = 2. 



Thus, for S 1 > 100000, 

If 1 I 0.30393 if «• ■ - J. 

'1/2' 



1 Z" 1 

(4.31) V- / 9w (uS/s)du < 

S <5 S 7l /2 



0.37273 if w = 2. 



For 5 < 100000, we proceed as above, but using the exact expression f|4.29[) 
instead of (EL"30j) . Note (T4T291) is of the form f s>m ,i(S) + f s , m ,2(S)/S, where both 
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fs,m,i(S) and fs,m,2(S) depend only on [S\ (and on s and m). Summing over 
m < S, we obtain a bound of the form 

1 f 1 

I g w {uS/s)du < G W (S) 



s<S 



1/2 



with 



G W (S) =K wA (\S\) + K Wj2 (\S\)/S, 
where K Wy \{n) and K w ^{n) can be computed explicitly for each integer n. (For 
example, G W (S) = 1 - 1/S for 1 < S < 2 and G W (S) = for S < 1.) 

It is easy to check numerically that this implies that (|4.3ip holds not just for 
S > 100000 but also for 40 < S < 100000 (if w = 1) or 16 < S < 100000 (if w = 
2). Using the fact that G W (S) is non-negative, we can compare G w (S)dS/ S 
with log(T + 1/N) for each T G [2,40] n -^Z (N a large integer) to show, again 
numerically, that 

(4.32) 



0.3698 log T ifu> 
0.37273 log T if w 



1, 
2. 



(We use iV = 100000 for w = 1; already iV = 10 gives us the answer above for 
u> = 2. Indeed, computations suggest the better bound 0.358 instead of 0.37273; 
we are committed to using 0.37273 because of ()4.3ip .) 
Multiplying by 6w /ir 2 a(w), we conclude that 



(4.33) Sx(U,W) 



if w = 1, 



(4.34) St{U,W) 



W 



W 



■Ho 



WU 



wu 



r 3/2 

5 - 08 « 3 / 2 ' W I 



if to = 2, 
(4.35) 



where 



0.22125 



Hence (by (JOS])) 
(4.36) 



if 1 < S < 40, 
if £ > 40, 



H 2 (s) 



0.15107 



if 1 < S < 16, 
if S > 16. 



0.22482 log T 
0.15107 log T 



if w 
if w 



1, 

2; 



moreover, £Ti(5) < 3/tt 2 , H 2 (S) < 2/vr 2 for all S. 

* * * 



Note. There is another way to obtain cancellation on /i, applicable when 
(x/W) > Uq (as is unfortunately never the case in our main application). For 
this alternative to be taken, one must either apply Cauchy-Schwarz on n rather 
than m (resulting in exponential sums over m) or lump together all m near each 
other and in the same congruence class modulo q before applying Cauchy-Schwarz 
on m (one can indeed do this if 5 is small). We could then write 

£ 2>(d) = - £ ^M^) = -E^)W+o(i)) 



m~W d\m 
m=r modq d>IJ 



m~W d\m 
m=r mod q ^ jj 



d<U 
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and obtain cancellation on d. If Uq > (x/W), however, the error term dominates. 



4.2. The sum S%: the large sieve, primes and tails. We must now bound 

2 



(4.37) 



S 2 (U',W',W)= 



U'<m<^ 



( log p)e( amp) 

W'<p<W 



for V = max(U,x/2W), W = max(S, W/2). 

From a modern perspective, this is clearly a case for a large sieve. It is also 
clear that we ought to try to apply a large sieve for sequences of prime support. 
What is subtler here is how to do things well for very large q (i.e., x/q small). 
This is in some sense a dual problem to that of q small, but it poses additional 
complications; for example, it is not obvious how to take advantage of prime 
support for very large q. 

As in type I, we avoid this entire issue by forbidding q large and then taking 
advantage of the error term 5/x in the approximation a = 2 -+- | . This is one 
of the main innovations here. Note this alternative method will allow us to take 
advantage of prime support. 

A key situation to study is that of frequencies on clustering around given 
rationals a/q while nevertheless keeping at a certain small distance from each 
other. 

Lemma 4.3. Let q > 1. Let a\, a 2 , ... ,ctk £ M/Z be of the form on = cn/q + Vi, 
< en < q, where the elements m £ M. all lie in an interval of length v > 0, and 



where at = a,j implies \Vi 
W,W >1,W> W/2, 



> v > 0. Assume v + v < 1/q. Then, for any 



(4.38) 



k 

E 

i=i 



E 

W'<p<W 



(logp)e(ay>) 



< min ( 1 



2q 



<j){q) \og{{q{v + v)) 



-l- 



(W 



W' + v- 1 ) ( lo §^ 

W'<p<W 



Proof. For any distinct i, j, the angles a^, ay are separated by at least v (if 



ai = a,j) or at least 1/q 



> 1/q — v > v (if Gtj a,j). Hence we can 



apply the large sieve (in the optimal N + 5 1 — 1 form due to Selberg |Sel91] and 
Montgomery- Vaughan |MV74j ) and obt ain the bound in (|4.38p with 1 instead of 
min(l, . . . ) immediately. 

We can also apply Montgomery 's inequality ([Mon68j, |Hux72j : see the expo- 
sitions in |Mon71j pp. 27-29] and [IK041 §7.4]). This gives us that the left side 
of (|4.38p is at most 
(4.39) 

/ \ 



-l 



E 

r<R 
\gcd(r,<j)=l 



{^{r)f 



J 



k 

E E E 

r<R a' modr i=l 
gcd(r,g)=l gcd(a'.r)=l 



E 

W'<p<W 



(logp)e((aj + a'/r)p) 



If we add all possible fractions of the form a'/r, r < R, gcd(r, q) = 1, to the 
fractions ai/q, we obtain fractions that are separated by at least 1/qR 2 . If v+v > 
1/qR 2 , then the resulting angles a« + a'/r are still separated by at least v. Thus 
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we can apply the large sieve to (14, 39ft ; setting R = 1/ \J [y + v )q, we see that we 
gain a factor of 
(4.40) 



> 



v (Mr)) 2 > HQ) (KrW > Hq) y 1 



logdqiv + v))- 1 ) 



ged(r,g)=l 

since Yld<R — l°g(-R) f° r all i? > 1 (integer or not). 



□ 



Let us first give a bound on sums of the type of S%(JJ,V,W) using prime 
support but not the error terms (or Lemma l4.3p . 

Lemma 4.4. Let W > 1, W > W/2. Let a = a/q + 0<(l/qQ), q<Q. Then 



(4.41) 



A <m<Ai 



(log p)e(amp) 

W'<p<W 

Ax-Aq 



< 



mm(q, \Q/2\) 

If q < W/2 and Q > 3.5W, the following bound also holds: 



{W-W' + 2q) Yj ( lQ gP) S 
W'<p<W 



(4.42) 



A <m<Ai 



(log p)e(amp) 

W'<p<W 
Ax-Aq 



< 



W 



log(W/2 9 ) 



W'<p<W 

If A\ — Aq < \5q\ and q < pQ, 5,p E [0, 1], the following bound also holds: 

2 

E E ( lo SP) e ( Qm P) 

(4.43) A <m<Ai lV'<p<W 

< (W - W + q/(l - Sp)) Yl ( l °SPf- 

W'<p<W 

The inequality ()4.42p can be stronger than ()4.42p only when q < W^/7.2638 . . . 
(if q is odd) or q < W/92.514 ... (if q is even). 

Proof. Let k = mm(q,\Q/2]) > \q/2]. We split (A ,Ax] into \(Ax - A )/k] 
blocks of at most k consecutive integers thq + 1, tjiq + 2, . . . . For m, m' in such a 
block, am and am' are separated by a distance of at least 

\{(a/q)(m - m')}\ - 0<(k/qQ) = 1/q - 0<(l/2q) > l/2q. 

By the large sieve 

2 



(4.44) Y 



a=l 



(log^)e(a(mo + a)p) 

W'<p<W 



<((W-W')+2q) ( lo grf 2 - 



We obtain (|4.4ip by summing over all \(Ax — Ao)/k~\ blocks. 

If Ax — Aq < \Sq\ and q < pQ, 5, p G [0, 1], we obtain (|4.43[) simply by applying 
the large sieve without splitting the interval Aq < m < Ax- 
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Let us now prove (|4.42j) . We will use Montgomery's inequality, followed by 
Montgomery and Vaughan's large sieve with weights. An angle a/q + a\/r\ is 
separated from other angles a! jq + a' 2 /r2 (r±,r2 < R, gcd(aj,rj) = 1) by at least 
1/qriR, rather than just 1/qR 2 . We will choose R so that qR 2 < Q; this implies 
1/Q < 1/qR 2 < 1/qnR. 

By Montgomery's inequality [IK04t Lemma 7.15], applied (for each 1 < a < q) 
to S(a) = ^2 n a n e(an) with a n = log(n)e(a(mo + a)n) if n is prime and a n = 
otherwise, 



1 

,(r) 



(4.45) 



} j (log p)e(a(m + a)p) 
W'<p<W 



< 



E 



a' modr 
gcd(a',r)=l 



^2 ( lo SP) e ( ( Q ( m o + a) + — ) p) 



W'<p<W 



for each square-free r < W'. We multiply both sides of (|4.45|) by (W/2 + 
(3/2)(l/qrR — l/Q)^ 1 )^ 1 and sum over all a = 0,1,..., q — 1 and all square- 
free r < R coprime to q; we will later make sure that R < W'. We obtain 
that 




1 

Q 



(logp)e(a(m + a)p) 

W'<p<W 



is at most 



(4.46) 




E E 

a=l a' modr 
gcd(a',r)=l 



1 

Q 



E 

W"<p<Vl/ 



(logp)e 



a (m + a) + — ) p 



We now apply the large sieve with weig hts [MV731 (1.6)], recalling that each 
angle a(mo + a) + a'/r is separated from the others by at least 1/rR — 1/Q; we 
obtain that f|4.46j) is at most J2w'<p<w(^°&P) 2 ■ (We- have followed a procedure 
analogous to that used in |MV73j to prove the Brun-Titchmarsh theorem.) 
Assume first that q < W/13.5. Set 



(4.47) 



R 



W 



(7- 



1/2 



where a = l/2e 2 -°- 25068 = 0.30285. 
R > 2. Moreover, for r < R, 



1 1 
Q ~ 3.5W 



< 



a 



1 



It is clear that qR 2 < Q, q < W' and 
cr/3.5 



a 



1 



< 



3.5 aW 3.5 qR 2 ~ qrR 
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Hence 



W 3 / 1 1Y 1 W 3 qrR W 3r W 
— + -(—--) <— + ^ =^r + ^n ^--2ct— 



2 2\qrR QJ ~ 2 2 1 - <r/3.5 2 2 (l - ^) i? 2 

W / 3<T rW\ W / 

t v 1 + w^itJ < t( 1 + ^ 



and so 



(-+9 

r<R 
gcd(r,(j)=l 



3/1 1 1 Kr) 2 



2 2 \qrR Q y 

gcd(r,g)=l 



V (1 + ri?- 1 )- 1 4rv > log + 0.25068; 



For R > 2, 

,-iwM r ) 5 

this is true for > 100 by [MV731 Lemma 8] and easily verifiable numerically 
for 2 < R < 100. (It suffices to verify this for R integer with r < R instead of 
r < R, as that is the worst case.) 
Now 

1 / W \ 1 W 

log R = - Hog — + log 2a) = - log — - 0.25068. 

Hence 

Tii + rR-r^y-^ 

^ y ' Mr) 2 & 2q 

r<R Yy ' H 

and the statement follows. 

Now consider the case q > VF/13.5. If q is even, then, in this range, inequality 
(|4.41|) is always better than (|4.42|) . and so we are done. Assume, then, that 
W/13.5 < q < W/2 and q is odd. We set R = 2; clearly qR 2 < W < Q and 
q < W/2 < W , and so this choice of R is valid. It remains to check that 

1 1 1 , W 

T H T > log • 

W,3fx_l_Y 1 W,3(±_iY 1 ~W & 2q 

2 2 \2q QJ 2 "T" 2 \4q QJ 

This follows because 

1 1 ^ , * 
r H r > log - 

1 4- 3 ( t _ JL}- 1 1 , 3 (t _ J_\~ 1 - 2 

2 "•" 2 \2 3.5/ 2 T 2U 3.5/ 

for all 2 < f < 13.5. 

□ 

The idea now (for large 5) is that, if S is not negligible, then, as m increases, 
am loops around the circle M/Z roughly repeats itself every q steps - but with 
a slight displacement. This displacement gives rise to a configuration to which 
Lemma 14.31 is applicable. 
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Proposition 4.5. Let x > W > I, W > W/2, U' > x/2W. Let Q > 3.5W. Let 

a = a/q + 5/x, gcd(a,q) = 1, \5/x\ < 1/qQ, q <Q. Let S 2 (U',W',W) be as in 

sm ). 

For q < pQ, where p £ [0, 1], 



(4. 48) S 2 (U',W',W) < (max(l,2p) (± + ±\+^ + 2qY £ Qogp) 

Ifq 
(4.49) 



lug pj 2 

W'<p<W 

Ifq < W/2, 



IfW> x/2q, the following bound also holds: 
(4.50) S2 (^< lW0 <(J + ^_) £ (W- 

If q < pQ, where p G [0, 1), and 5^0, 
(4.51) 

2q/cj ) {q) 



S 2 {U',W',W) < min 1, 



I0g(^(g+2W) 1 



( X W X X \ v-^ . 2 

KM + T + 4Tw)q + 5^7) EN ■ 



W'<p<W 

The trivial bound would be in the order of 

S 2 (U',W',W) = (x/2logx) Yl ( l °gP) 2 - 

W'<p<W 

In practice, (|4.50p gets applied when W > x/q. 

Proof. Let us first prove statements (14.490 and (|4.48j) . which do not involve 8. 
Assume first q < W/2. Then, by (Q2l) with A = U', A x = x/W, 

V 9 J <K?) log(W/2?) w ,^ w 

Clearly {x/W - U')W < (x/2W) ■ W = x/2. Thus (IQ9l holds. 

Assume now that q < pQ. Apply (|4.4ip with yl = U', Ai = x/W. Then 

^'• W '- W) £ ( ,. m lTl,Prt-) + ') (W ~ W ' + 29) E (l0grt2 ' 
Now 



+ lV - W' + 2q) 



q ■ min(l, (2p) 1 ) 

fx \ ~W — W' / x \ 

< ( 777 - U ') -——— + 2mzx(l,2p) [—-U') + W/2 + 2q 

~ \W I gmin(l,(2p)- 1 ) v ' HJ \W J ' H 

gmm(l, (2/9j 1 ) Vv 

This implies t[4~gg]) . 

If W > x/2q, apply i4l3| with 5 = x/2Wq, p = 1. This yields (ITOjl . 
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Lastly, assume 5/0 and q < pQ. Let Q' = x/\5q\; then Q' > Q. For any mi, 
m2 with |mi — m®] < Q' — q, 

| (mi — 7742) • <V X + S^/ 35 ! — Q'$/ x — 

We split the range U' < m < x/W into |~(ac/W — U)/(Q' — q)~\ intervals of 
length < Q' — q; in each interval, the conditions of Lemma 14.31 are fulfilled with 
v = (x/2W) ■ 5/x and v = \5q\/x. We obtain a total upper bound of 

S 2 (U',W',W) < min (l, 2,1 



<j)(q) log^ + f))- 1 ), 



W'<p<W 

Here W - W + v~ l = W - W + xjqb < W/2 + x/q5 and 



Moreover, 



W + x^\ ( 1+ x/W-U\(W +( j\ r x/2W 



2 qSJ \ Q'-q J ~ V 2 /V (1 - 

< -7T + Q + 



2 ^ 4(1 -p)Q' 2W(l-p) 
x W x x 



\Sq\ 2 4(1 -p)Q 2(1 -p)W 

□ 



5. Totals 

Let x be given. We will choose U, V, W later; assume from the start that 
U < x/2e, 2 • 10 6 < V < x/4 and UV < x. Starting in section [5T2| we will also 
assume that x > x = (2e) 2/5 • 10 25 = 6.891459 . . . • 10 26 . 

Let a £ R/Z be given. We choose an approximation a = a/q + 5/x, gcd(a, q) = 
1, q < Q, \d/ x \ < 1/qQ- We assume Q > \fex and Q — max(4£/, x/U). Let S^i, 
5*7,21 Sq be as in (|2.16p . with the smoothing function 77 = t]q as in (jl.4p . 

The term So is because V < x/4 and 770 is supported on [—1/4, 1]. Consider 
w = 1 and w = 2. 



5.1. Contributions of different types. 



5.1.1. Type / terms: Si 1. The term S^i can be handled directly by Lemma 13.41 
with po = 4. (Condition (|3.1ip is valid thanks to (|2.4p .) We will always have 
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Q/A> U . Thus the contribution of Sj i is bounded by (|3.12p and ()3. 13p : 
(5.1) 

Sii = xrtq)<Kw/ q +) y Mm) lQg ^ 

q w/q + ^-^ m m 

* ' m<M /q- 

gcd(m,wq~)=l 

+ ( ^l^) ^ { _ §) _ y Mkgd ^ . y /x(ro) 

W / q ft d ) m<M 0/q - m 

gcd(m,wq~ )=1 



, i • Q , , Uw ( ' lQcex ex , , 
+ 0<[ - log 2e<? log x H log — — log — - + 1 • 8 log 2 



7T 7T \ ?7ii; 

+0s ((-r + " (i " ) -s- log — J ■ U - s?) ' 48 J ■ 

where we use (|2,3p . and where we set Mo = min(J7, |_a;/(|&/|)J /•«;), <7 + = gcd{q,w), 
q- = q/ q + and c = 1 + 1/(4 log 2). Note that M = U (because \5q\/x < 1/Q 
and Q/2w > U) and 



fi(d) log d JO if w = 1 or q is even, 

d | — if w = 2 or q is odd. 

a|u>/q+ v ^ 



Recall (Ramare; ([£gD ) that 

(5.2) V ^l g^I = 0<(1.4 + ^-M^-)). 

z — ' m m 

m<!//5 
gcd(m,«!5 - )=l 

Let 

r?i(i)=log f-^t^(t) 
By (H3D and the fact that x/U > 4, 

s 04(^(1,^)). log 1 



Now 



77! (i) = log -r] + log -777 % = log -J/o - log 777^-%- 



Hence the first two lines of (|5.ip equal 

Xfi(q) 4>(w/q + )^ s-^ Km) U/q~ 
t—Vo{-o)- > log 

q w/q + m m 

m<U/q~ 
gcd(m,wq~ )=1 

[ xfijq) ( (f)(w/q+) y fi(d)log(d) \ y fi(m) 

q \ wlq + ^ d J ^ m 

\ 1 d\w/q+ / m<U/q- 

gcd(rn,wq~)=l 
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By (|2.5p . (|2.7p and (j5.2j) . the absolute value of this is bounded by 



43 



(5-4) - min ( 1, ( ^ + TT^T + ( bg ^ ] A 



if u; = 1 or g is even, and by 



2x 



if Tif = 2 and g is odd, where 



, 2q/H J ] if w(q-) 2 < U/e 2 ^l 
(5.6) Ai= { l ° s ^F7 ~ 

1 otherwise. 



Since (a) log 2 < 1.4 and (b) u>g /<ft(wq x ) = q/4>{q) if w = 1 or g is even, (c) 
wq~ /4>{wq _1 ) = 2q/4>(q) if to = 2 and g is odd, (|5.4p and (|5.5p are both bounded 

by 



We obtain 



/ 59.1792a; ex \ 
c 1 ) + 1.7651(7 log 2eg log x + 1.7651C/u; I log — '— log — h 1 

v ' \ Uw Uw J 

+ 3.5683 aw log — — + — ^ , 

\ xq U e I 



where Ai is as in (|5.6p . 



5.1.2. Type I terms: Sj^- The case q < Q/V. If q < Q/V, then, for v < V, 



va n / d \ va _ / 1 \ 
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and so va/q is a valid approximation to da. We can thus estimate Sj 2 by applying 

Lemma 13.31 to each inner double sum in Sj 2 (12.161) . We obtain 

(5.9) 

S/,2 = A(t>) 53 fJ-(m) e((va) ■ mn)rj(mn/(x/v)) 

v<V m<U n 

gcd(«,tu)=l gcd(m,u))=l gcd(n,m)=l 

xfj,(q v ) 4>(w/q+) v-^ /i(m) 



E A (< 



v<V 
gcd(u,tu)=l 



m<M v /q v 
gcd(m,ui(ji7)=l 



mu 



+ E A ^- Q < ( , ^! 9 : )fTH -^-^i^u 



gcd(u,io)=l 



xq v /vw 



8 2vr 2 / 



/ a,, , [/tt> / . , , 2ecx , 16ecx \ 

A(v)|V|i- ( -log2e 9 „ + I 2(1 + e) log-— + lop 1 

v<V ^ n \ v 



Uvw J 



gcd(t!,llj) = l 



53 A(b)(2Hi + |i/|i).( 



, / 3wcx/v 3(1 + e) x/w C/io 



gcd(i>,iu)=l 



Q/v + 2e Q/v log M„ 



where we write g„ = qj gcd(w, <?), g+ = gcd(g„, w), q v = q v /q+, c = 2 \vh/W \l + 1 
and My = max(min(C/, Q/2wv), 1). We bound the first sum as in (|3.17p . using 
(|2.7P ; we obtain 



^— ^ ' DO.. nn/n+ ^ ' 



V<V 
gcd(ti,«))=l 



y: a( 

gcd(l),M)) = l 

<^^A 2 l g(^) 



:cd(g,v) 



m<M/q v 
gcd(m,wq~)=l 



III 



V 



v w/qt 



log 



by (pTTUj) . where 
(5.10) 



■^#2L if t«(?- ) 2 < My /e 2 ^(^) , 
1 otherwise. 



The second sum is at most 
(5.11) 



1 1\ , . U „ 
-- w ja(w)w-\r, |i 



53 A(u)(f7 + q)v + 53 A(t>)u • £/gcd(?;,<7) 
gcd(«,g)^l 
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By {USD, 

Mv){U + q)v < 1.03884— (U + q). 

v<V 

Proceeding a little coarsely, 



£ A(v)vgcd(v,q)<J2 E 0ogp)-p Q -p min(Q ' 

v<V p\q a:p a <V 

gcd(u,g)=l 

< ^(logp) .pLlogV/log,j p « p (,) < 2Vq\ogq. 
p\q 

Thus (|5.1ip is at most 

a( W ) W ( 1.8535 ^ + V + 7.137^ log^ . 
\ gx x J 

We can use (|2.12p and partial summation to bound 

E. , . /„.., Xl 2ecx . 16ecx\ 
A( W ) (2(1 + C ) log _ + log — J 



v<V 



2e 2 cx , 16e 2 cx 



< 1.03883V 2(1 + e) log — — + log 



UV UVw 

We bound other error terms in (|5.9j) using (|2.1ip . bounding q v from above by q 
when needed. 

We obtain 
(5-12) 

IS | < ^Mz^a 2 i og( y g) + ( 2(1 + 6) log ^ + log 

<7 \ (71/ l/ku; 

+ a(wW ( 1.8535 ^ + g V 2 + 7.137^-^ log g 

\ qx x 

+ 1.7658* log 2eq + 30.8119-^- + 11.3223^-^ log 

Q/V Qt/V Q 

where My = max(min(C7, Q/2wqV), 1). We will choose e > later. 

The case q > Q/V. Lemma 13.51 gives us (|3.14j) and (|3.15j) . which we use 
directly. 

5.1.3. Type II terms. The quantity to bound here is (|4.4p . that is, 

[■X/U JTIA fx/U JTIA 

(5.13) 4^ v^XU W) ■ S 2 (U, V,W)— + 4j ^ Sl (U,W)-S 3 (W) — , 

where Si, S2 and S3 are as in (|4.5p and (|4.6p . We bounded Si in (|4.33p and 
(Q41 . S 2 in Prop. g3]and S3 in g^I). 

We first recall our estimate for Si. In the whole range \V,x/U] for W, we know 
from (|4T33l) and (P4"|) that Si(U,W) is at most 



3 x 1 x/WU 

if u; = 1, 



, r1r , 2 X Wo/nx3 X /X/W 
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if w = 2, where 

'5.08 iiw = l, 



1.27 if w = 2. 



We have better estimates for the constant in front in some parts of the range; 
in what is usually the main part, (|4.33p and (|4.34j) give us constants of 0.22125 
and 0.15107 instead of 3/V 2 and 2/vr 2 . Note 5.08C(3/2) 3 = 90.5670... and 
1.27C(3/2) 3 = 22.6417. . . . We should choose U, V so that the first term domi- 
nates. For the while being, assume only 

(5.16) f/^5.10 5 -^; 
then (|5.14p and (|5.15|) give 

(5.17) S^U,W) < i^^, 
where 

A + 9 /- 5 f 1 < 0-4321 ]£w = l, 

4 + 4S«8 < 0.2347 ifw = 2, 

and this will suffice for our cruder estimates. 

The second integral in (|5,13[) is now easy to bound. By (|4.6|) . 

S 3 (W) < 1.0171a; + 2.0341W < 1.0172a:, 

since W < x/U < x/10 6 . Hence 

rx/U dW f x/U I — — dW 

'-'if • 1M72x w 



r 1 / dw t 

4 1 ywv) ¥ <4/ v 



v 



where 

'5.3035 if«j = l, 
3.9086 if «; = 2. 



Kl0,9 = 8 • y/1.0172 • K^i < 



Let us now examine 5 , 2, which was bounded in Prop. 14.51 Recall W' = 
max(V, W/2), U' = max(*7, x/2W). Since W > W/2 and W > V > 117, we can 
always bound 

(5.18) Oogp) 2 < ^(logW). 

W'<p<W 

by (EH]). 

Bounding S2 for 5 arbitrary. We set 

Wo = min(max(20(/, V),x/U), 

where 9 > e is a parameter that will be set later. 
For F < W < W , we use the bound (|Q5|> : 

S 2 (tf',W,W) < (max(l,2p) f|- + ^ + ^ + 2g) • ±W(log W) 

A \ /W \ W 2 logW 
< max -,/0 — + 1 )x log W H -2 h qW log W, 



4q 



where p = q/Q. 
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If W > V, the contribution of the terms with V < W < W Q to (|5.13[) is (by 
I5.17P bounded by 



Wo 
V 



/ x f Po (W A , xxr W 2 logW \ 



(5-19) +Kw>3] J— + qx J v —^dW 



where /?o = max(l, 2p), 



+ gx( (log W f /2 - (logy 



«w,2 = 4^2*^,1 < 

and 

« m ,3 = 4:y/^w,l < I 



I 3.7183 if w = l, 

1 2.7403 if u; = 2 

2.62921 if io = 1, 

1.93768 if io = 2. 



We now examine the terms with W > Wq. (If 6q > cc/J7, then Wo = U/x, the 
contribution of the case is nil, and the computations below can be ignored.) 
We use (14491) : 

s,v,w,w) < (^j-^ + ^j-^) . Iww. 

By \/a + 6 < -y/a + Vb, we can take out the q/(f)(q) ■ W/log(W/2q) term and 
estimate its contribution on its own; it is at most 

r-x/U 



[ / — JL lw2 l ° gW — 
J Wo Y^W' 4>{q)' * \ogW/2q W 



Hi) J Wo VW \ V log W/2q J 



Now 

r /u i / io S 2o , r IWq i 

/ -==* -, hl—dW < j2q\og2q / dt. 

J Wo y/W\logW/2q J ma x(e,V/2q) y/t\ogt 

We bound this last integral somewhat crudely: for T > e, 

f T i / T 

(5.21) / dt < 2.34 



y/Uogi ~ ' V logT 
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(by numerical work for e <T <Tq and by comparison of derivatives for T > Tg , 
where T = e^" 2 / 2 - 3 )" 1 = 2135.94 . . . ). Since > e, this gives us that 



x ' u 1 / / log2g \ nir 



and so (|5.20p is at most 



2^/^(1 + 1.15/ l0g2(? 



<j){q) y " y \ogx/2Uq y y^" 
We are left with what will usually be the main term, viz., 



'■■■'''/ i >• logW \ dW 



( 5 - 22 ) 4 \ Si(U,W)- \w 



W 



A(/){q)logW/2qJ W' 



which, by (|4.33p and (|4.34|) . is at most xj y (j)(q) times the integral of 



1 

W' 



, jx/WU\ logW 



•\ piWi ' W ' 4 V 17 y log^/2g 
for W going from Wq to x/U, where H w is as in ()4.35p and 



362.2684 if io = 1, 
90.5671 if w = 2. 



= 4k Wi0 C(3/2) 3 < 

By the arithmetic/geometric mean inequality, the integrand is at most 1/W times 



B + B- 1 -4H w (x/WU) B~ l x/WU 8 log2g 
(5.23) h —t-k w m — — h -- 



U 2 log W/2q 

for any B > 0. We will choose /3 later. 

The first summand in (15. 23p gives what we can think of as the main or worst 
term in the whole paper; let us compute it first. The integral is 

r*/u p + p-i. AH w (x/WU) dW f x / UWo 6 + /3" 1 ■ 4:H w (s) ds 



I 

Jw 



(5.24) 



v v 2 + 2/3 ) io & UW 11 W 

by (|4.36p . where 




„■„ - W Jx 2 s 

§_ , *™,6 \ l X :r _ -i 

1 ^ 2/3 I 1U & C/W W ~~ ' 



Thus the main term is simply 



0.89928 if w = 1, 
0.60428 if w = 2. 



(5.25) ('^ + ^6 > \_^ = iog X 



2 2/3 y y^y ° [W 

The integral of the second summand is at most 



2 U h W*/ 2 WA V U 



P '~ FT I ~77PTF) — P ■ K w,4y - 
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By (|5.16p . this is at most 

/r 1 V2 ■ ht 3 • k wA < /r 1 ^, 

where 

V2k w4 Jo.5124 if u; = l, 
KwJ ~ 1000 - jo.1281 if to = 2. 

Thus the contribution of the second summand is at most 

a— 1 ^ 

P ^7 



The integral of the third summand is 

P f x/u log2g dW 

2j Wo T^gw/^W- 

If V < 20q < x/U, this is 

P f x/u log2g dW B f x/2Uq 1 eft 

- log 2q 



(5.26) 



2 J 20q logW/2q W 2 b J e logt t 

P ( x 

= g lo S 2< 7 • (l°g l°g 2^ ~ log log e 

If 20(7 > x/U, the integral is over an empty range and its contribution is hence 0. 
If 20q < V, (^26]) is 

P f x / U \og2q dW P\og2q f x / 2U l 1 dt 



\ogW/2q W 2 J v/2q logt t 

p\og2q 



(log log — - log log V/2q) 

3\og2q ^ , \ogx/UV\ 
log 



(Of course, log(l + {\ogx/UV)/{\ogV/2q)) < Qogx/UV)/QogV/2q); this is 
smaller than (logx/tft^/log 2g when V/2g > 2q.) 
The total bound for (|5.22p is thus 

(5 - 27> vk ' (" ' G log ^ + *) + r 1 log + K - 7 ) ) • 

where 

'log2g (log log ^ - log log ej if F/20 < g < x/{29U). 



(5.28) $ 

Choosing P optimally, we obtain that (|5.22p is at most 



log^iiog^i + ^i^g-) if«<v/29. 



(5 ' 29) Tfe v v log ^ + *) v""' 6 log m + 2K " J 

where $ is as in f|5.28|) . 

Bounding S2 for \S\ > 4. Let us see how much a non-zero 5 can help us. It 
makes sense to apply (|4.51 j) only when \6\ > 4; otherwise (|4.49p is almost certainly 
better. Now, by definition, \5\/x < 1/qQ, and so \5\ > 4 can happen only when 
q < x/2Q. 
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With this in mind, let us apply (|4.5ip . Note first that 

yi > 2/\5g\ > mm(Q,W) 



x 

W\ 



+ 



2WJ 



Q ^ W 



\Sq\ 



Recalling also (j5. 18[) . we see that ()4.51j) gives us 

2q/<f>(q) 



S 2 (U',W',W) < min 1, 



log ™ap 



( 



x W x x \ l TTr/1 TTT . 



where p = q/Q. Note also that \6q\ < x/Q. 

The contribution (to (15.13j) ) of the terms W/2 and x/(4(l — p)Q) can be set 
aside (by y/a + b < yfa + v6) and bounded somewhat crudely; it is at most 



x/U 



W + 



2(l-p)Q 



\ ^ max (l,™ log 



mi n(Q,HQ 



WlogW 



dW 
~W 



(5.30) < 2^/K~i 2^= + j=l = log ■ X 



VU y/2(l -p)Q b UV 



max 

v<w<§ 



logW 



max [l^Miog—^. 



min(Q,iy) 



Since Q > x/U, mm(Q,W) = W for all W < x/U. Since t (->• (logt)/(logt/c) is 
decreasing for t > c, it follows that 



max 
V<W<-, 



logW 



min(Q,WQ 



< 



/ 2g log!/ 



if |<5<?| < V (as is the necessarily the case if Q > ex/V). However, if \Sq\ > V/9, 
we use the bound 



max 

y<W<jj 



\ogW 



max I 1, 2— lQ g — [^i — 



< 



^log^N 



min ( 1, ^^jp log I 



The contribution of x/2(l — p)W can also be set aside; since in it the values 
of W close to V dominate, we choose to bound the contribution by 



x/u / x ( x \ Vi w dw 



V 



4^ 



w , 1 



2y/T 



P Jv 



x/U 



TIogTF 



, 2-v/ Ktu 1 

yj K w,l 



dW 



v 



< 4 



iogT + v/i/iog v) x, 
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where we use the estimate 

3 \^ogW 1 f°° Vlogn + logl/ 



,V^gV , 1 , / 2 



What is left to estimate then is the main term, coming from x/qS. Similarly 
to before, we define Wq = max(V, 0|<5g|), where 6 > e will be set later. We bound 
the main term by 2x/ yj $4>{q) times the integral of 



<x/WU\ logW 

Hu,4~ 



from Wo to x/U, plus (if Wq > W) the integral 

r 6 \ &q \ / x 1 dW 

(5.32) ^ Sl(U ,W). w -W^W w . 

This contribution of (|5.31|) is bounded (much as in (|5.24|) ) by 



lo § 77TI7- + $ K ->6 log + 2 KwJ ) , 



where 



(5.33) $ 



log |^| log l + fef^X ifN<W 



iogV/1%, 

log \5 q \ ( log log ^ - log log e ) if v/e < \5 q \ < x/eu. 



We have Wq > W if and only if 9\5q\ > V; in that case, we must also consider 
the integral (|5.32[) . which we bound as follows: 



f e \ s i\ / x 1 dW , x f 6 ^ sfioe. W 

A l f^-w\-2 wlogW w^^^ml ^ dW 



6 vW V 



It is time to collect all type II terms. Let us start with the case of general 5. 
We will set 6 > e later. If q < V/29, then \Su\ is at most 
(5-34) 

x If x ( iogx/t/y\ 



x , / \ogxl\JV\\ / x 

log — + log 2q log \1 + j j (^ i6 log — + 2*,,, 



+ 2k Wi2 , -T7-T 1 + L15 1/ i TTTTT" "7= + «w,9 



^ y logx/2UqJ ' VV' 
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If V/29 <q< x/29U, then \S n \ is at most 
(5.35) 



log 2g \ x 



+ 2kw,2\ rm\ 1 + L15 W i 777FT -7= + («tu,2\Aog 2 #<7 + K w,9 



\ogx/2Uq y/U ' v.^v~*-, ■ 



+ ^|((log2^) 3 / 2 -(logF) 3 / 2 )^ 



<?x, 



where we use the fact that Q > x/£7 (implying that po = max(l, 2q/Q) equals 1 
for q < x/2U). Finally, if q > x/26U, 



\Sn\ < (K W}2 \/2]ogx/U + K Wt9 )—= + K w $y/\ogx/U 



v vu 

(5.36) / \ 

+ ^((iogx/c/) 3 / 2 - (bgy) 3 / 2 ) [-^= + V^y 

Now let us examine the alternative bounds for |<5| > 4. If \5q\ < V/9, then 
|<Sjj| is at most 
(5-37) 

2x /A x 1 ip m A logx/[/F\\ / . a; 

V V log uv + log 1591 log I 1 + bgWj J ^ log uv + 2Kw > 7 



4^/^Tt 



X 



+ V m \Vu + 78(Q^) log uv) ' V log f/n ' " 

If |£g| > then |5//| is at most 
(5.38) 

2x 



l0g E7»j% + l0g M bg logg J V K - 6 bg c7^ + 2K - 7 J 



+ ^^((iog^|^|) 3 / 2 -(io g y) 3 / 2 

6 ^/\6q\ V 



x/bgi 7 + yr/ioiy) + 



^,2 ( J_ 1 1 i g I ^log 6>|<5g| 

V2 { vu vm^) oguv )' rrj, im 



min [l,±?f-\ogO 

where p = q/Q- (Note \S\ < x/Qq implies p < x/AQ 2 , and so p will be very small 
and Q — q will be very close to Q.) 

5.2. Adjusting parameters. Calculations. We must bound the exponential 
sum ^2 n A(n)e(an)r](n/x). By (|2.15p . it is enough to sum the bounds (|5.8p for 
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d5H2D for 5 7j2 (or (13+5]) if g > Q/V), (|532]) - (|08D for 5//. We will now 
see how it will be best to set U, V and other parameters. 
Usually, the largest terms will be 

/ . . , 2e 2 cx 16e 2 cx 
(5.39) 1.8337L7Vuj 2(1 + e) log — — + log 



UV ° UVw 

(from (|5,12p . type I) and 
(5.40) 



x 



log -J- + log \25 iq \ log ( 1 + *° Su y ] ] (W 6 log -J- + 2k Wi7 



l/v 



(from (fCTD and (IOTP , type II; here S = max(l, |<5|/4), <5i = 1 if \S\ < 4 and 
<5i = \S\/2 if |<5| > 4). For now, we will think of e as being (say) 1/4. 
We must choose UV. We simpljH set UV = x/t, where 

t = Ky/qEl, 

and K will be optimised for the range of greatest concern (q ~ 2 • 10 6 , 
4 • 10 5 ). What we must minimize is x times 

j (23.07416 - 1.8337 log u; + 6.41795 log t) 



+ 



log t + log \28 iq \ log [ 1 + - lQg ^ J J (K wfi log t + 2k wJ ) 



\/<M<?)<5o \ \ V log 



\28 iq \ 



For 6 = Si = 1, q = 2 • 10 6 , cp(q) ~ 4 • 10 5 , and V = 10 15 /^ (a good value to 
examine, as later arguments will show), this is (to five significant figures) 

1 (69.63208 + 6.41795 log K) 

K-V2TW 

1 I " 1 7.25433 + log K 

+ 4 / log if + 7.25433 + 15.20180 log 1 + 



V4 • 10 5 V ' V log 13.12236 

• \J 0.89925 log K + 7.24806 
if w = 1, 

'136.72211 + 12.8359 log K) 



K • V2 • 10 6 



1 / / 7.25433 + log K 

+ W log K + 7.25433 + 15.20180 log 1 + 



V4 • 10 5 V ' V log 13.12236 



• v/0.60428 log A" + 4.56485 

if w = 2. (Larger choices of w are also possible, but preliminary computations 
make it doubtful that they give better results.) 
The choice 

[26 if«j = l, 



K = K„ 



81 if w 



seems best or nearly best. We settle for K\ = 25, K2 = 49.8 to give ourselves 
some leg-room. 



*^An alternative method would be to approximate log(l + (log x/UV)/(\og V/2q)) by 
(log x/UV)/(\og V/2q) and then use the non-principal branch of the Lambert function. 
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We set UV = x/t = x/K w ^q5o, at least for q considerably smaller than V . 
(The exact value of q at which the choice of t will change will be set later.) Then 
([Q9]) plus (HlHD is at most 



x 



X 



3.21888 + log |25ig| I ^ + log I 1 + 1 



(5.41) vW)\ V 2 V log 



v 



\25 iq \ 



■ a/0.44963 log qSi + 3.61917 + -=(0.11919 log q8 x + 1.74931) 



for w = 1, and 



3.91202 + log |2<5iq| - + log 1 + 



(5.42) vTWfaJ\ ~ ' ! V 2 V log 



v 



• a/0.30214 log g + 2.54516 + —^==(0.11919 log g<5i + 1.87642) 

for u; = 2. We will work with w = 2 from now on, since it is now clear that it is 
the better choice: compare the terms under the second root in (I5.4ip and (15.421) . 

Now we must decide how to choose U, V and Q, given our choice of UV. 

Our tactic at this point will be to choose Q ~ x/Ay (say) and use the Sip 
estimates for q < Q/V to treat all a of the form a = a/q + 0<(l/qQ), q < y. 
The remaining a then get treated with the (coarser) Sip estimate for q > Q/V, 
with Q reset to a lower value (call HQ'). If a was not treated in the first go 
(so that it must be dealt with the coarser estimate) then a = a'/q' + S'/x, where 
either q' > y or S'q' > x/Q = Ay. The value of Q' is set to be smaller than Q 
both because this is helpful (it diminishes error terms that would be large for 
large q) and because this is now harmless (since we are no longer assuming that 
q < Q/V). 

5.2.1. First choice of parameters; q < y. The largest items affected by our choices 
at this point are 

x 2UVw 2 

(5.43) 56.6115^-^ log+ — - — , 1.7658? log 2eq (from S If2 , e = 1/4), 



In addition, we have a relatively mild but important dependence on V in the 
main term (|5.40p . We must also respect the condition q < Q/V, the lower bound 
on U given by (|5.16p and the assumptions made at the beginning of section [5] 
(e.g. Q > x/U). It is from this point onwards that we assume x > xq = 3.1 • 10 28 ; 
this will simplify many of our choices. We also set w = 2, and so K w = 50. 

The following choices then make sense: we will work with q < y, where 
(5.45) 

x 1 



y = 2x - 3 , Q = — = -x - 7 , x/UV = K w Jq~6\ < 50j2y~= lOoVaP*, 
4y 8 

X 0A5 jjy x 4x o.55 x oa 
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where, as before, 

Si 



1 if \S\ < 4, 
|(J|/2 if |5| > 4. 



(Note that 5/x < 1/qQ implies \q5\ < x/Q = 4y, and thus q5\ < 2y.) 
Since Q = x/4y, 

Q V = 7777 = = K w — ^/qd l = — ——yJqOiy 

UV x 4y 4 Ay 

= > J° = 1.1048 ...q>q. 

4 4 • 2 3 / 2 V H y ~ 32V2 

Moreover, U/(x/UV) > x a3 /400 > 8.81 • 10 5 > 5 • 10 5 , and so (f57l6l) holds. 

It is easy to check that 
(5.46) 

U < x/4.33, 2 • 10 6 < V < x/4, UV < x, Q > ^/ex, Q > max(J7, x/U), 

as stated at the beginning of section ([5]). We will also find it convenient that 

V V x2 4 /25 x o.i 

> > ' = -2— = 2.599 . . . > 1, 

2eq ~ 2ey ~ 4ex[j- 3 lOOe 

' 5 ' V V t - 1 

V>V= ^Q_ = 1.2996... >1. 
e|«5g| ~ 4ey 200e 

The first type I bound is 

+ 1.7651 • y log(2ey) log x + 0.88255 • x a45 (log(118.4x a55 ) log(2ex a55 ) + l) 
+ 3.5683 (^.^(V/V-)^) 

< _L_ . f 2.4 + 21 °/( 2 ^ 55 ) , I + 0.073*°«, 
■Him \ max (2, log 



where (52 = max(l, (7r/2y3)*) and Ai is as in ()5.6|) . It is easy to show that, for q 
fixed, 

(5.48) 21 ° S < 2l °' 55) 



. _ , T 0.45 

max 2, log =5-3- 



attains its maximum at x = (8e 2 q 2 ) 1 ^ 0A5 , and this maximum equals 
log ((8eyf 55 /°- 45 /2) = f log ? + (f + flog 2) ; 

of course, (15.480 is also at most log(2x 0,55 ) in general. Thus 

x /22 11 \ 

\Sia\< — * -min — log q + 6.693, — log x + 3.094 + 0.073x' 

(f>(q)°2 V 9 20 / 
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The second type I bound is (by (j5. 12[l . with e = 0.06) 

x ( 48 \ x 

|S/, 2 | < - • min 1, — — X 2 log(Vq) + -= • (1.68373 + 0.11443 log (gfc)) 
q V (2vr5)V y/qJx 

/ fl /IK / (\ AK. \ 

0.45/4 

+ 7.137- 1 



+ 6 1.8535 



4 I 4 



4^,0.55 \ 2 ^0.45 , 

L ' 1 - - ^log(2x - 3 ) 
50x/^ SV ' 



50V^] 



,0.85 



32x 

+ 1.7658 • 2x a3 log (Aex - 3 ) + 30.812 

25y qd\ 



84g 0.85 g / 50 

+ 11.3223 — ^ log — Q^y- 



< 



<t>{q)8\ 



7 A 2 log(Vg) + 0.0045^- + O.OOlx 06 
i- q z 6\ 



' 75^^ S x°- 7 /8 



x / , , r N 38.406 log x + 71.042 

^= • 1.63303 + 0.11076 log (qS{\ + 

qdi V x 

where A 2 is as in (|5.10p . We simplify the last line by 



, , . 38.406 log x + 71.042 1 m ^ 
I ^ < 1-63303 + ^ < 1.771. 

x o 



1 «™ ^ 38.406 log x + 71.042 
1.63303 H ^rrs 

Note now that 



m ( ■ (tt Q \ A Q K ™ U 

My = max I mm I U, — — ] . 1 ] — 



4V I6y yyul 16-8 



50 



)\X 



,0.15 



Thus 



(5.49) 



0(9) 



\ 1 M/ W 21 °g ^ 
A 2 log(Fg) < ? — 

9 max M2, log 

Z1 ° g l 50V7T 



2 log 



4x a55 v /q 



50 



max (2, log ^#)- max (2, log 

For x < (e 2 • 128g 3 / 2 /25) 1/0 - 15 , this is at most 

/ 4j^yg \ ^ 4(e 2 -128g 3 / 2 /25)°- 55 /°- 15 ^ 

log — < log — < 61ogo + 10.796. 

& V 50 7 ~ & ^ 50 j - 

For x > (e 2 • 128g 3 / 2 /25) 1 / - 15 , the expression (|5.49p is of the form c(t + a)/(t + b) 
(with t = logx), and hence is monotonic on t (and thus on x); hence it adopts 
its maximum either at x = (e 2 • 128g 3 / 2 /25) 1 / - 15 or at x = oo. Therefore, its 
maximum is 



max 



0.55 
0T5 



,log 



/ 4 / eM|g^ \°- 55 / - 15 



V 



50 




= 6 log q + 10.796. 



Hence 



|S7,2| < —TTE2 min ( lo § 



<j>(q)8 



4x°- 55 ^q 



50 



, 6 log q+ 10.796 



+ 



-|= • (1.771 + 0.11076 log(g<Ji)) + 0.0045— ^ + O.OOlx 06 . 
/qdi q z di 
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For type II, we have to consider two cases: (a) \5\ general, and (b) |<5| > 4. We 
have q < V/2e and \5q\ < V/e thanks to (]5.47p . We apply (]5.34p with 9 = e, and 
obtain 



\Sn\ < ' 



3.91202 + 



■ ^0.30214 log q + 2.6202 

+ 5.4806, / [1 + 1.15./ lQg .^ rr | + 3.9086— 



< 



25 

X 



^3.91202 + C x , q log 2q- x/0.30214 log q + 2.6202 



+ 4.8971, —^x - 8 + 19.543X ' 8 , 
V 0W) 



where we define 



C x ,t := - + log | 1 + ; -g^g- 



2 V log 



25t 3 / 2 



for t such that 25i 3 / 2 < x°' 55 /2e (as is always the case here, thanks to q < V/2e). 
In case (b), we use ()5.37p (again with 9 = e) and obtain 



2x i 

S u \ < 7=== • a/ 3 - 91202 + Cx,5 19 log N • v/0.302141og 5 lQ + 2.6202 

4^/Ku,,! /- — /— — \ X 



+ v ^1/logF) + 



+ «»*^\^=+ , 1 log(100x - 15 ) • J^p- • x 

' \<l>(q)\VU y/8(l - y/Q)Q J ^ogV/Ay 

< X ■ ^3.91202 + C xAq ■ log 2faq- \/0.302141og 5^ + 2.6202 
\A)<A(<?) 



Now note the fact f |RS62l Thm. 15]) that g/0(g) < F{q), where 

2.50637 



(5.50) F (q) = e 7 log log q + 



log log q ' 



It is easy to show that F(q) < F (y) = F (2x a3 ) and that x H> x om /y/g (2x - 3 ) is 
increasing for x > xq (indeed for x > 1). Hence 



max j 19.5 .13 + 1 .897 1 A j 77.339 + 4.5493W ) < 3.33575x a0fl 
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Hence, the total bound for q < y is at most 

x 



\A)</>(<?) 

(5.51) 



• ^3.91202 + C XtSiq • log 25^ ■ ^0.30214 log 5^ + 2.6202 

+ ■ (1.771 + 0.11076 log (qSi)) 

x (11 1 .76 

+ , , . r9 min — log x + - log q + 0.5683, — log q + 17.489 
0(g)o| \10 2 9 



'2 

Sit 



0.0045-Ar + 3.3358X ' 85 , 



where 



Cr 



X.t 




log 50^ ^ 

± \ __n kk I 



AOg 25*3/2 / 



(Jo = max(l, |<5|/4), 



62 = max I 1, — —8 



[l if|<5|<4 
1 1<5|/2 if |*| > 4, 



2^3 

Note that i 1— )■ Cj, j is an increasing function of t, and thus attains its maximum on 
[0, 2y] at t = 2y. (We recall that 5±q < 2y.) Both x h-> and x i-)- C Xj ( (fixed 

i) are decreasing functions. For x = x = 3.1 • 10 28 , we have C^^-io 6 = 1-196 . . . , 
C x>y = 2.463 ... and C x>2y = 3.021 .... 

5.2.2. Second choice of parameters. If, with the original choice of parameters, we 
obtained q > y = 2x 0,3 , we now reset our parameters (Q, U and V). We will 
do so in such a way as to minimize the terms of size about UV. (This makes 
sense because q is large, and thus the terms proportional to are already fairly 
small.) It would not be sensible to make UV much smaller than Q, since the error 
terms proportional to Q would then dominate. 
We choose 

V = 2-x°-\ C/ = 500x - 35 , g = £ = * 

U 500 

Thus (|5.16p holds (as an equality). The conditions in (|5.46[) also hold. 

Write a = a/q + S/x for the new approximation; we must have either q > y 
or \5\ > ^y/q, since otherwise a/q would already be a valid approximation under 
the first choice of parameters. 

Let Ai be as in (|5.6H . Clearly Ai < 1. Hence, by fl5.8B . 



1 - i_4-W^1 



( 59 1792x ex 
+ 1.7651Qlog 2eQlogx + 2- 1.7651C/ flog— — — log — 



+ 3.5683 




< 



Aj- ( -|-r + 0.65 log x - 5.5 I +2.961x 2/3 . 
m \<t>{q) ' ) 
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4UV / 2UV , 8ecx 



By (IXT51) with e = 0.001, 



|5/, 2 | < -2log(^) + 41og2 log log. 

8.008C/T/ / ^ 2ecx , 

+ 4 log 2 log 2C/y log — — + log 2ecx 

IT \ UV 

+ 4 log 2 • ( — logff/V + Q) log 16cx + — log 2eQ log Q 
\ ir ir 

+ ( ^_ 4) ^l(,„3 88 3f + !^.(5) 2 + ^ (logQ)Q 



/12cs 3003cx, 8tfV\ 
1 log 

V Q Q Q 



+ 4log2 • ( + — - — log^^J log2J7y 



< —= (0.3 log x + log q + log 2) + 20.563x L 



Now we must estimate 5//. There are two cases: (a) \6\ < 4 and (b) \S\ > 4. As 
we said before, we must have q > y or \6q\ > Ay, and thus in case (a) we certainly 
have q > y. Since y = 2x 0,3 = V, this implies, obviously, that q > V > V/2e. 
Hence we can treat (a) by applying (|5.35 j) with 8 = e, provided that q < x/(2eU). 
Let us bound some of the terms in (|5.35|) . First of all, 
(5.52) 



2K w ,2«/^-r I 1 + 1-15, 



log 2q 



<f>{q) y y logx/2Uq J y/U 

< 2k w 2 \ F (^tt) ( 1 + 1.15 4 /log 4? J ; < 0.04421x - 85 , 

W ' Z ]/ \2eUJ\ V 6 e[/y V500x - 35 



where F (q) is as in (I5.50P . (We are using the fact that second expression in (|5.52|) 
divided by x - 85 is decreasing for x > 1.5 • 10 14 .) In much the same way, we obtain 
that 



,5,3) ^ 3 (f^ + <y**w- <y*vr* \ * < , 4088I .. 85 



using the easily checked fact that the quotient of the left side of (|5.53p by the right 
side of (|5.53p is a decreasing function of x for x > 2.4 • 10 28 . Since t i— > (logi) 3 / 2 
is convex-down for t > 1.7 and t h-> i/e'/ 2 has its maximum at t = 2, 
(5.54) 

3 n m i/2 lo g 2e 9 - log^ 



((log2eg) 3 / 2 -(logy) 3 / 2 )^<|(logy) 1 / 



3 1/a log e 2 V- logV / 21ogy 

< -(log K ) 7 . X = 3\ — — — X. 

- 2 K ^ ' ^eV/2 V eV 
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Hence 
(5.55) 

|S "' - Tfe ' v /log *k + log 21 log log 4 ' v /o - 60428 log m + °- 2562 



+ (j^^logV + K w<2 \/logx/U + K^ )9 ^ 

+ 0.04421x - 85 + 0.4088X 085 



< 



y^y v log ^ + log29loglog m ~ V - 60428108 2^ " a34808 



+ ( 1.7898^/1^+3.21681 - 

V log x J 



where we are using the fact that 1 1— > v£ is convex-down. 

Let us now consider the other subcase of (a), namely, q > x/{2eU). Then we 
use (|5.36p . and get that 



(5.56) 



\Sll\< {K Wt2 sf\ogx/U + K W! g/V2)x°- 85 + K w>3 \fiogX~/U 



U 



+ 



K w ,3 



((logx/Uf/i-ilogVf/^t^ 



exU x 
+ 



U 



< (2.7403 y/\ogx/U + 0.09122 Vlog a; + 2.8653)x 

( i 4.2366 \ n « 

< 2.3006Vlog2; + 2.86 53 - - x 85 . 
V Vlogx/ 



0.85 



We now study case (b), namely, |<5| > 4. We know that \5q\ > Ay; we also know 
that \8\/x < 1/qQ, and so q < x/{\5\Q) < x/AQ = U/A and \Sq\ < x/Q = U. We 
treat (b) by applying f|5.38j> : the second line of (|5.38|) can be treated just as in 
(j534"|) . We obtain 

(5-57) 

2x / x x I x 
\S n \ < — 7=== A /log 777773- + log \Sq\ log log 77777 • J/t W|6 log + 2k^ i7 



v1¥(g)V e^N 



C/|<5g| 



eC/|<5 9 | 



. K »,2 /logy 



4Q 



+ K w> 2 



<P(q) 



u 



8Q ( 1 - m 



aW + -7== ) + ^,9 1 

VlogV^ 
\ 

^=^= log 



ylog eU ■ x 



J 



< 



2x 



log 777777 + 1^1 1°S 



U\5q\ 



Kwfi log 



eC/|<5g| 



+ 2k 



w,7 



+ 1.95801 ydog 2x - 3 + 3.1522 + 



1.37038 
v/log 2x° 



,0.85 



.3 
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Now we must add all the contributions. We aim at a bound of the form 
C(log log x log x)x ' 85 . The terms with xjqb\ in front give a total of 
(5.58) 



0.95 log x + log q + -f- + (log 2 - 5.5) < - • 1.54024 < 0.77012x 



x 



„0.7 



qb\ V. °~ ' ™ * ' 0(g) ' v ~ & ~ '"V y 
Hence 

| + \Si i2 \ < 0.77012x 7 + 2.951x 2/3 + 20.563x 085 
< 0.074922(logloga;)(logx)x a85 . 



The term in (|5.55p with x/ \J (p(q) in front can be bounded by 



r( q ) L Q . , Q . Q 9 

X V ~q~V g 2eq + log2gloglog Km, > 6 g 2^ + ™' 7 ' 

this is decreasing for q E [y,Q/2e], and thus is bounded by its value at q = y, 
which is at most 

0.160884(log loga;)(log2;)x a85 . 



The other terms in (15. 55ft add up to at most 

0.061521(log logx)(log3;)x a85 . 

We can easily bound all of ()5.56|) by 

0.076427(log logx)(log3;)x a85 . 

Recalling that \5q\ > Ay in case (b), we see that the expression in (|5.57p with 
2x I y / \5\(j)(q) in front is at most 

0.157099(log logx)(logx)x a85 , 
and the other terms in (|5.57p add up to at most 

0.044793(log logx)(log2;)x a85 . 
Hence (|5.55|) is worst, and the total bound is at most 
(5.59) 0.297327(loglogx)(log:r)x a85 , 
where we use the fact that 0.297327 = 0.074922 + 0.160884 + 0.061521. 

5.3. Conclusion. 

Proof of main theorem. We have shown that \S v (a, x)\ is at most (|5,5ip for q < 
2x 0,3 and at most (|5.59|) for q > 2x ' 3 . It remains to simplify (|5.51|) slightly. Let 
g = 2-10 6 . 

= 3.91202 + q EO , qo - log 2 g0 
9 0.30214 log go + 2.6202 

By the geometric mean/arithmetic mean inequality, 



^3.91202 + C xAq • log 25^ • ^0.30214 log 5^ + 2.6202 
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is at most 



\ ( ^(3.91202 + C X)Siq ■ log 2hq) + Vp ■ (0.302141og 5 x q + 2.6202) 

< 3.2423 + i ^C^, + VP" 0.30214^) log 2^ 

< 3.2423 + | 0.4091 + 0.28151 log | 1 + /° g 5 °^^ | | log 25 x q. 



log 



25(5 ig )3/ 2 



Hence, for q < 2x°" 



._ , 3.2423 + log 2£ig 1.6943 + 0.1108 log 2^? 
5,(a,a;) < x + = x 

/ll 1 76 \ x 

+ min — log x + - log g + 0.5683, — log q + 17.489 — — = 
V10 2 9 J (j)(q)5% 

+ 0.0045-A7 + 3.336x 85 , 
Siq 2 



where 



R x , t = 0.4091 + 0.28151 log f 1 + ]° g 5 ° ^ 



log +i 



25t 3 / 2 



□ 



In practice, the simplification of (|5.5ip leads to a loss of less than 0.05 percent 
in the worst case (for x ~ xq), and much less than that in the sensitive range 
near q = 2 • 10 6 . 



6. The integral over the minor arcs 

We have obtained bounds for S^x, a), where a is in the union m C S 1 of all 
minor arcs. We must now use these bounds to estimate the total contribution of 
the minor arcs 



(6.1) / Srj (x, a) (S Vl (x, a)) e(—Na) da 

J m 

to the integral over all of S , which gives us the (weighted) number of represen- 
tations of N as a sum of three primes: 

J ^S V0 (x,a)(S m (x,a)) 2 e(-Na) da = ^ ry (J^j rjl (^) m (^7) • 

" ni+n2+ns=N 

Here 770 is as in f j 1 . 4 [) . i.e., it is the smoothing we have been using all along; r/i is 
an arbitrary smoothing. (It would make sense to choose rji so as to maximise the 
major-arcs contribution, or so as to make that contribution easy to compute.) 
One would usually set x equal to N times a constant to be optimised. 
The traditional procedure here is to bound (|6.ip by an ' £2 estimate: 

S„ (x, a) (S Vl (x, a)) 2 e(— Na) da < (max5„ (x, a)) ■ \S Vl (x, 

„ ' (■ •. m 
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This introduces a factor of |£„.(A,a)|| ~ ClogiV, C a constant. We will use a 
technique due to Heath-Brownjj based, in turn, on a lemma of Montgomery's we 
had the chance to use before. Asymptotically, this replaces Clog N by a constant; 
it is this fact that ensures that only a bounded number of L-functions will have 
to be checked for the major-arc treatment. For the sensitive range consisting of 
values of x not much larger than xq = 3.1 • 10 28 , the improvement will be more 
slight. 

Lemma 6.1 (Heath-Brown). Let S(a) = Y^ n a n e(an), where a n 6 C has compact 
support contained in (\/x, oo) and a n = whenever n is divisible by any prime 
p < r. Let 



3.2) 



m 



U U 

q<r a mod q 
gcd(a,g)=l 



a r a r 
q yq q yq 



for r > 1, y £ (2r 2 ,2xr 2 ]. Then 



-i 



3.3) 



|5(a)| 2 da < 



E 



/ 



There is a minor difference between this and the formulation in [Taol Lemma 
4.6] - namely, here we allow the length of the intervals to vary. Compare to the 
idea of the weighted large sieve in [MV73] , which, like the method here, takes 
advantage of the unequal spacing of Farey fractions. 



Proof. Apply Montgomery's lemma jlK04i Lem. 7.15] to the sequence a n e{an) 
with a£R/Z, u(p) = l. We obtain 



<P(q) 



eian, 



a mod q 
gcd(a,g)=l 



^a n e((a + a/q)n) 



for any product q of distinct primes p < \fx. Letting R = y/y/2/r < \fx and 
summing over all square- free q < R coprime to all primes < r, we obtain 



(6.4) 



a n e(an) 



< 



z E E 

q<R a mod q 
q sq-free gcd(a,ij)=l 
p|g=>p>r 



y^a n e((a + a/g)n) 



Where L = Eq<R: P \q^p>r Kq?/Hq)- 

Let <7i,<72 < r, ai < qi, a 2 < q2 with gcd(ai,qi) = 1 and a\/qi / a 2 /q2- For 
any q[,q' 2 < R, a[ < q[, a' 2 < q' 2 with gcd(a-,g-) = 1 and gcd(<?-, <?ig 2 ) = 1, 
the angles a\/q\ + a\jq\ and 02/^2 + o! 2 /q' 2 are separated by at least \jq\q 2 R 2 . 



^Personal communication (2006). The technique has since been described in [Taol §4]. 
Granville and Balog previously suggested a plausible alternative, namely, using mollifiers as 
in |Bal90] , 
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Since R = ^Jy/2/r, we know that \lq\q2R 2 = 2r 2 jyq\q2 > (2/y) max(r/gi, rjq^). 
Hence the intervals 

ai a\_ t_ ai a£ _r_\ / 02 r a 2 | gg | r 

91 ?i yqi Qi q[ VQiJ' W q' 2 yqi Q2 g' 2 vqi 

are disjoint. This implies that all translations m^^ + a/q (with q < R square-free 
and coprime to all primes < r, and a coprime to q) are disjoint. 
Integrating and summing (|6.4p over all intervals in m y r , we obtain 

/ \S(a)\ 2 da<\ EE/ a \S(a)\ 2 da 

Jm v> T q<R a mod? Jm v,r+^ 

q sq-free gcd(a,q)=l 
p\q=>p>r 



<j[ |5(a)| 2 da<lE 



a n \ 2 . 



□ 

Later we will need explicit lower bounds on the sum on q in (|6.3p . By (|4.40p 
(an estimation standard in sieve theory; see, e.g., (ircOll pp. 63-64]), 



e %>-uH)-y:} 



q<R p<r v r/ d<R 

p\q=^p>r 



Also, 



w n f 1 - 1) ■ e 



1 log # + 7- 
P <r v ^ £^ e^logr + j^ 

for J? > 1, where we are bounding the denominator by [RS62, Thm. 8]. Recall 
R = y/y/2/r in our case, and so logi? = (log(y/2))/2 — logr. 

If R < r, then we obviously cannot do better than the trivial bound 

E ^ 

q<R 
p\q=i-p>r 

If R > r, we can also use the bound 

m(<?) 2 1 log-R 1 

6.8 V ^->1+ > ->l + l0g-^ rj, 

z — ' 0(q) P logr (logrr 

q<R r<p<R V & V 6 7 

p|g=>p>r 

(for r > 286, by [ RS62 , Thm. 5]), which is better than (|6.5p when r < R < r 4 or 
so. Further refinements are possible if R > r 2 ; (|6.5p is optimal only asymptotically 
as (log i?)/(log r) — >• oo. 

Lemma 6.2. Let f, g : {a, a + 1, . . . , b} — > Wq , where a, b £ Z + . Assume that, 
for all x G [a, 6], 

(6.9) E /(")<^). 

a<n<x 
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where F : [a, b] — > E is continuous, piecewise differentiable and non- decreasing. 
Then 

b rb 

/(n) • g(n) < (maxg(n)) • F(a) + / (max^(n)) • F'(u)du. 

^— ' n>a I ri>u 

n=a ~ Ja - 

Proof. Let S(n) = J2m=a f( m )- Then, by partial summation, 

6 b-l 
(6.10) ■ 9( n ) < S ( b )9(b) + £ S(n){g{n) - g{n + 1)). 



Let h(x) = max x < n <b g(n). Then h is non-increasing. Hence (|6.9p and (|6.10p 
imply that 

b b 

<£/(«)&(«) 

n=a n=a 

6-1 

< S(6)/i(6) + £ S(n)(h(n) - h(n + 1)) 



n=a 
6-1 



< F(b)h{b) + £ F(n)(h(n) - h(n + 1)). 

n=a 

In general, for a ra G C, ^(x) = X^a<n<x a « ana ^ ^ continuous and piecewise 
differentiable on [a,x], 

rx 

a n i ? (x) = A(x)F(x) — / A{u)F' {u)du. {Abel summation) 



a<n<x 



Applying this with a n = h(n) — h(n+l) and A(x) = Yl a <n<x a « = h{a) — h{[x\ + 
1), we obtain 

6-1 

J2F{n)(h{n) - h(n + l)) 



6-1 



(%) - h{b))F{b - 1) - / (/»(a) - /i(M + l))F'(u)du 



,•6-1 

h{a)F{a) - h{b)F{b - 1) + / /»([uj + l)F'(u)du 

J a 

h(a)F{a) - h{b)F{b - 1) + / h(u)F'(u)du 

J a 

h(a)F(a) - h(b)F(b) + [ h(u)F'{u)du, 



since h([u\ + 1) = h(u) for u ^ Z. Hence 

6 



^/(n) 5 (n) = h(a)F(a) + [ h(u)F'(u)di 

„—„ J a 



□ 



Proposition 6.3. Let 5(a) = ^ n a n e(an), a n G C, {a n } in L 1 (Z). and m y)r 
6e as in Lemma \6.1l Let T : R/Z — > C. Lei ro,ri G Z + , ri > rp > 1 and 
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y 6 (2r 2 , 2xr 2 ]. Assume that a n = whenever n is divisible by any prime p < r\. 
Let g : [ro, r\] — > M + be a non-increasing function such that 

(6.11) max \T(a)\<g(r) 

ae(R/Z)\m y , r 

for all r £ M, ro < r < r%. Let H : [ro,ri] — > M + be a non- decreasing, continuous 
function, differentiable almost everywhere, such that 

(6.12) 1 I \S{a)\ 2 da < H{r) 

2-i \ a n\ Jm yt r + i 

for all r G [ro,ri], and, moreover, H{r\) = 1. 
TTien 

' ; \T{a)\\S(a)\ 2 da 

< g{r Q ) ■ (H(r ) - h) + T g{r)H'(r)dr, 



En Kl 



■'•() 



where 

E„ l» 



(6.13) /o = ^ ^ , 2 / |S(a)| 2 c2a. 



Proof. For ro < r < r%, let 



/(r) = =-_ / 15(a)) 



2 



TOy,r+l\niH,' 



and 



En \ a n 



/(^l) = V I 12 / 

Z^n \ a n\ J(«./Z)\m y , ri 



Then, by (IBTTT]) . 

?4l5 / |S(a)| 2 |T(a)|da < £ /(r)*(r). 

By (EH, 



^ /(r) = =J-^ / \S(a)\ 2 da 
, <r<x l^n l a ™l ■/m»,*+i\ 

m y, r o 

\S(a)\ 2 da ) -J < H(x) - I 



rn<r<x 

(6.14) , 

/ 1 



EnKI 



for x G [ro,ri). Moreover, 

£ /(r) = [ \S(a)\ 2 

,<vo, 2^nl a n| J(R/Z)\ my , ro 

= ( V I | 2 / |5(a)| 2 ) -/ = l-/o = ^(ri)-/ 



ro<r<ri 
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We let F(x) = H(x) — Iq and apply Lemma 1631 with a = r , b = r±. We obtain 
that 

f(r)g(r) < (max g(r))F(r ) + / (maxg(r))F' (u) du 

*■ — * r>ro l rn r>u 

r=ro u 

< g(r )(H(r ) - J ) + g(u)H'(u) du. 

□ 



Lemma 6.4. Let x > x = 3.1 • 10 28 . Lei 

(6.15) R x , r = co + ci log ( 1 + 1 , 

where Co = 0.4091, ci = 0.28151. T/ien t/ie function 

is decreasing on [l,2x - 3 ] provided that f3 > 0.175. 
Proof. Write 

log 50^ , , s \/r 

/i( r ) = " -oITi h{r) = 

log ^372 log4rVloglogi 

We need to show that 

+ log (l + /i(r)))' f[{r) 



c + ci log (1 + h(r)) (1 + /i(r))(c + c x log(l + f x (r))) 

is less than 

m /? 

h{r) r 

for r G [l,2x - 3 ]. We can write /i(r) = (loga)/(31og 6/a), where a = 2500r and 
6 = (5000x - 55 ) 2 / 3 . Therefore, 

df i . . d log a log 6 11 log b 

(r) = 2500- -2- r = 2500- 



dr da 3(log 6 — log a) 3(log 6 — log a) 2 a r3(log6/a) 2 

and so 



d/((r) . 1 (Iog6)/(Iog6/a) 



2 



(1 + /i(r))(c + ci log(l + /x(r))) 3r / loga \ A / loga 



log 6 11 



2r ^ J + | log 6) ( t + log (l + ggfc) ) 2r »(«)0°8 & ) ' 
where a = (log a) / (log b) and 

w(a) = (1 - a)((l - a) + 1/2) U/ Cl + log (l + ° ^ . 

The assumption x > xq implies that (log a)/ (log b) < 0.949. It is easy to check 
that q — >• w(a) decreases on [0, 1) and that w(0.949) > 0.0963. Thus, 

w(a) log b > 0.0963 • 29.732 > 2.863. 

Hence 

11 1/3 

< < 



2rw(a)(logb) 5.726r r 
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for Q > 0.175 > 1/5.726. □ 

Lemma 6.5. Let x > xq = 3.1 • 10 28 . Then the function g : (e, 2x 0,3 ] — > M + given 

by the following is non-increasing for r > 43/ 

(6.16) 

g(r) = ((3.2423 + R x , r log 2r) \J F{r) + 0.1108 log 2r + 1.6943) 

0.85 



IT 

/ll 1 76 \ Fir) 

+ min — log x + - log r + 0.5728, — - log r + 17.4945 J —^x + 3.336x 
V 10 2 9 

where R X)r is as in 16.15\) and 

, 2.50637 
(6.17) F(r) = e 7 log log r + 



log log r 

Proof. Taking derivatives, we can easily see that 

ill 76 \ log log r 

min I — log x + - log r + 0.5728, — log r + 17.4945 ' 



10 ° 2 ° '9° IT 

is a decreasing function for r > 10, 

N 1 

r -> 0.1108 log 2r + 1.6943)-= 
is a decreasing function for r > 1, 



Vlog log r 
r -> 3.2423 v 6 _ 
\/r 

is a decreasing function for r > 14. These statements remain true if log log r is 
replaced by F(r), since F(r)/ log log r is a decreasing function for r > 1. 
We can also check easily that 

Vlog log r 



log 2r- 



r 0.325 



decreases for r > 43. Together with Lemma \6A\ this implies that 

r ^ R lor 2r V1 ° g bg r - lor 2r V1 ° g bg r 

r ita. )t . log zr ^_ - r0 .i75 log /r r 0.325 

is a decreasing function for r > 43, and so we are done. □ 

Lemma 6.6. Let x > x = 3.1 • 10 28 . Let 
(6.18) 

g 2 (r) = ((3.2423 + R X)2r log4r) \J F{r) + 0.1108 log 4r + 1.6943^ 



11 , 1 . 76 



\ F(r) , 



„0.85 



+ min — log x + - log r + 0.5728, — log r + 17.4945 -^^x + 3.336x 
\10 2 9 / r 

where R X)T is as in t6.15\) and F(r) is as in \6.11^ . Then, for r > 35000, 

<? 2 (2.25r) < g(r), 

where g is as in \6.1b}) . 

Proof. By Lemma HT41 R x ,4.5r < 4.5 - 175 i2 Ijr . Since r > 35000 and F(r)/loglogr 
is a decreasing function, 

log 9r v / F(2.25r) log9rVloglog2.25r V^25 

/ — < — ; = — < 1.1528 < , i . 

log2r A /r(r) log2rVloglogr 4.5 - 175 
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This implies that 

(R x , 2ri log4n)^5^ < (i? Xir log2r)^^ 

for ri = 2.25r; the rest is simpler. □ 
Lemma 6.7. For any q > 1 and any r > max(3, q), 

" <F{r) 



for F as in {6.11^ . 

Proof. Since F(r) is increasing for r > 27, the statement follows immediately for 
q > 27 by [RS621 Thm. 15]: 

< F{q) < F(r). 



For r < 27, it is clear that q/4>(q) < 2 • 3/(1 • 2) =3; it is also to see that F{r) > 3 
for all r > e. □ 

Theorem 6.8 (Total of minor arcs) . Letx > xo, xo = 3.1-10 28 . Let S^a^x) be as 
in $1.1}) : let rjo be as in 1 1-4\ ), and let r/i : [0, oo) — > [0, oo) be a bounded, piecewise 
continuous function of fast enough decay (namely, lim^ 00 (\/i log t) ■ r]f(t) = 
and | log(t)r]i(t/x)\i < co). Assume that the series in hl.l\) converges absolutely 
for rj = T)\. Let m^ ,. be as in Lemma \6.1\ Let rg > 16000. 

Ifr < (x/18) 1 / 4 , then 
(6-19) 

\Srj (a, x)\\S m (a, x)\ 2 da 

(R/Z)\m y , ro 

< J(g{ro)-{H{r + l)-L )+ J\(r)H'{r + l)dr^j • ^ A(n) 2 »ft(n/x) 2 



+ max(5(ro), 0.2974 (log log x)(log x)x 



0.85^1/2 



• J 0.7131 ^ • m 2 (j^j dt + 1.03883|r/i|^ • n logn. 

where g(r) is as in A6.16\) . n = (x/18) 1 / 4 , Lq is as in A6.13]) . 

'« 7 k)+i) i 



(6.20) H(r) = min 

and x(r) = (x/18) 1 / 2 /r. 
For ro arbitrary, 

1 



<(l-/o) 



logx(r) + 7 - ^ ' l + Jog i2£i£M L 



log r (logr) 2 



U / s / / M2 / 5 »» (a, a?) \\S m (a,x)\ da 

21) f 

max (5(r ),0.2974(loglogx)(logx)x ' 85 ) i/r < 2x a3 , 

0.2974(log log x) (log x)x - 85 if r > 2x a3 . 
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Proof. Let Q = x 0,7 /8, as in the Main Theorem, and y = x/9. Let a G (M/Z) \ 
m yr , where r > tq. There exists an approximation a = a/q + 5/x with q < Q, 
\5\/x<l/qQ. 

There are three cases: 

(1) q < r. Since a is not in m^, then, by definition (|6,2p . \5\/x > r/yq = 
9r/xq; that is to say, \5\ > 9r/q, and, in particular, \5\ > 9 > 4. Thus, 
for (5q, <5i and 62 as in (jl.3p . 



$oq = -jq, Siq = -jQi W?) = ^=|%- 

Recall that q/4>(q) < F(5oq) by Lemma [6.71 we note also that 5oq < 
(x/qQ)/4 ■ q = x/4Q = 2x ' 3 . Therefore, the main theorem implies that 
\S V0 (x,a)\ < g2($oq), where 52 is as in f)6. 18|) . Now 5oq = \5\q/4: > 
{9r/q)q/A = 2.25r > 2.25tq. By Lemma 16.61 (applicable since 2.25ro > 
2.25 • 16000 > 35000) and Lemma EB 



52(<5o<?) < 9 (J^) ^ 9{r). 

Therefore, 

(6.22) \S V0 (x,a)\<g(r). 

(2) r < q < 2x ' 3 . Then, by the main theorem and Lemma |6.5| 

(6.23) \S V0 (x,a)\<g(q)<g(r). 

(3) q > 2x - 3 . By the main theorem, 

(6.24) \S V0 (x,a)\ < 0.2974(loglogx)(logx)x a85 . 

We set r = ro and bound \S V0 (x, a)\ for a G m^^ as above. In consequence, 

^ 77-T2 7 1 \2 I \S V0 {a,x)\\S m {a,x)\ 2 da 

<^ — rr~^ — 7 — r~v? ( / \S Vl {a,x)\ 2 da\ ■ max \S Vo (a,x)\ 

J2n A (n) 2 Vl(n/xy \J im) \ mytro l Vn J a G (R/Z)\m w , ro V ° 

<(l-Jo)- max \S Vo (a,x)\ 



< (1 - Io) 



max (sr(ro),0.2974(loglogx)(logx)x°- 85 ) if r < 2x 03 , 
0.2974(loglogx)(logx)x - 85 if r > 2x a3 . 



Assume now that ro < (x/8) 1 / 4 . Let r\ = (x/8) 1//4 . Write 
m (a,x) — ^ (logp)e(ap)7?i(p/x), 

p>ri 

m (a,x)— A(n)e(an)r/i(n/x) + y^ A(n)e(an)?]i(n/x). 

n non-prime n<r\ 
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By the triangle inequality (with weights \S r)0 (a,x)\), 



\Sr! (a, x)\\S m (a, x)\ 2 da 
-2 J / \S m {a,x)\\S j)711 {a,x)\ 2 da. 

j=1 V J(TSL/Z)\my,r 



Clearly, 



\S V0 (a,x)\\S 2 , Vl (a,x)\ 2 da 

X V,rQ 



< max \S„„(a, x)| 1//2 • , / / \So m («, x)\ 2 da 
ae(W/z)\m y , ro V Jr/z ' V 



< max \S TK Aa,x)\ 1 ^ 2 ■ / A(n) 2 r]i(n/x) 2 + 'S^ A(n) 2 r]i(n/x) 2 . 

ae(R/Z)\m y ,r \ ■ ~^ 

u y n non-prime n<ri 

Let = sup w>z r]i(w). By [RS62, Thm. 12], [RS62, Thm. 13] and partial 

summation, 

y~] k{n) 2 r] 1 (n/x) 2 < ^ A(n) 2 r/T(n/x) 2 

n non-prime n non-prime 

/ \ 



< 



n) 2 



(rjl 2 (t/x))' dt 



n<t 

\n non-prime / 



/oo _ 
(logt) • 1.4262 y/i (rji 2 (t/x))' dt 

< o.7i3i r^.wiu 



^2 A(n) 2 r]i(n/x) 2 < ( lo g r i) X] A ( n ) 

n<r'i n<r'i 

< 1.03883|t/i|^ • nlogn. 



It remains to bound 

)ll^l 

We must check that the conditions of Prop. 16.31 hold for S(a) = Si tVl (a,x), 
T(a) = S vo (a,x). It is clear that y = x/9 E (2r 2 ,2xrQ]. We defined Si,^ in such 
a way as to let it be supported on integers coprime to all primes p < r\. By 
Lemma 16.51 g{r) is non-increasing on [ro,ri]. In particular, g(r) > g{r\). It is 
easy to show that 

g (n) = g ((x/8) 1/4 ) > 0.2974(logloga;)(logx)x a85 

for all x > x . Hence (l6T22|) . (I6T23|) and (^21) imply that 

\S m (x,a)\ < g(r) 

for all a E (K/Z) \ m^.. 
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Inequality (|6,12p (with H(r + 1) instead of H(r)) holds by Lemma 16.11 and 
([63 |) -(|6T5 ]) (with R = ^/y~/2/r). It is clear that H : [r + l,n + 1] M+ 
is non-decreasing, continuous and differentiable almost everywhere, since it is 
the minimum of three functions satisfying these three conditions and crossing at 
finitely many places. It remains to check that H{r\ + 1) = 1. This follows easily 
from ([(HOj) and the fact that >c{r\ + 1) < r x + 1. 

We apply Prop. 16.31 an d obtain 

\S V0 (a,x)\ ■ \S 1;Vl (a, x)\ 2 da 

/Z)\m v ,r 

< (g(r )-(H(r + l)-I ) + H g(r)H'(r + l)dr) ■ A(n) 2 m (n/x) 2 . 

□ 

Let us discuss briefly how to read the statement of Theorem 16.81 The last 
two lines of (|6.19p are negligible. If ro < r\, the leading term in H'{r) on the 
range [ro,ri] (save for the rightmost part, where H'(r) = 0) is proportional to 
1/rlogr. For g(r) as in (|6.16p . this implies that J" 1 g(r)H'(r + l)dr is bounded 
by (essentially) a constant times g(ro)/ log ro- In contrast, the main term g(ro) ■ 
H{r$ + 1) is (in practice) a constant times g{ro). 

The size of the negative, and hence helpful, term —g{ro)Io depends on the 
quality of major-arc estimates, but it is always within — C<7(ro)(logro)/logx and 
for a constant C. 

For x = x = 3.1 • 10 28 , ti6^0D gives us that H{2 ■ 10 6 ) = 0.85578... and 
H(A ■ 10 6 ) = 0.92115 .... Numerical integration gives 

rn f 0.02316... forr = 2-10 6 , 

0.01773... forri = 4-10 6 , 



g(r ) • H(r + 1) + H g(r)H'(r + l)dr 

J ro 



which is only slightly better than the bounds given more naively by (|6.2ip . namely, 
g(r ) = 0.02379 ... (for r = 2 • 10 6 ) and g(r ) = 0.01782 ... (for r = 4 • 10 6 ). 
The eventual major-arc help —g(ro)I would be the same in each case. 

The true power of the improvement of (|6.20p over (|6.19j) becomes clear only 
for x much larger than xo- As x increases and r remains fixed, H(r) asymptotes 
to 2e 7 (logr)/(logx). The denominator cancels the factor of logx implicit in 
J2 n A(ra) 2 77i(n/x) 2 (see (|6.19p ). It is because of this that it will be enough to 
check major arcs contributions up to a constant ro, even as x — > oo. See Table 
[21 which gives values of 

(6.25) ^g(r ) • #(r + 1) + J g(r)H'(r + l)drj ■ log x 

for variable x and ro = 4 • 10 6 . (Prolonging the table after 10 80 is unnecessary, 
since H(tq) is already equal to the first term within min in (I6.20p . and so <?(ro) • 
i^(ro + 1) • logx cannot grow.) 

There is a middle range (roughly: x between r 6 and r 10 ) for which Theorem 
16.81 can be refined further by the inclusion of a fourth term in (|6.20p . The reason 
why the middle term in (|6.20p is sometimes better than the first term is that 1 + 
Yl r <p<z l lv is a better approximation to J2 n <z:p\n^ P >r 1 / n than ( lo g z )l ( e7 lo § r ) 
when z is between r and r 4 (or a little further). When z is between r 2 and r 4 , we 
can also gain something by considering terms of the form I/P1P2, v < Pi,P2 < z/r, 
P\P2 < Z- This can be used to improve slightly on the middle range in Table [2j 
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3.1 • 

10 29 
10 30 
10 32 
10 34 
10 36 

1Q 40 
10 50 
1Q 60 

10 70 
1Q 80 



Bound on minor-arc integral 



llF" 1.163699. 

1.162418. 
1.153979. 
1.124607.. 
1.090788.. 
1.059498. 
1.012654.. 
0.975229 . 
0.997130.. 
1.042510.. 
1.028066., 

Table 2. Size of the bound (1(1251) for r = 4 • 10 6 



Appendix A. Computational work on n = p\ + p2 + P3 for n small 

In order to obtain an no such that every odd number n < uq is the sum of 
three primes, we need only combined known estimates on the binary Goldbach 
problem (n = p\ + P2) and on prime numbers in short intervals. Estimates on 
primes in short intervals are generally based on numerical verifications of the 
Riemann hypothesis up to a certain height, i.e., checking that all (non-trivial) 
zeroes z of the Riemann zeta function up to a height H (meaning: < H) lie 
on the critical line = 1/2. 

The height up to which the Riemann hypothesis has actually been fully verified 
is not a matter on which there is unanimity. The strongest claim in the literature 
is in [GD04] . which states that the first 10 13 zeroes of the Riemann zeta function 
lie on the critical line dt(z) = 1/2. This corresponds to checking the Riemann 
hypothesis up to height H = 2.44599T0 12 . It is unclear whether this computation 
was or could be easily made rigorous; as pointed out in [SDlOt p. 2398], it has 
not been replicated yet. 

Before [GD04], the strongest results were those of the ZetaGrid distributed 
computing project led by S. Wedeniwski |Wed05j ; the method followed in it 
was more traditional, and should allow rigorous verification involving interval 
arithmetic. Unfortunately, the results were never formally published. The state- 
ment that the ZetaGrid project verified the first 9 • 10 11 zeroes (corresponding to 
H = 2.419 • 10 11 ) is often quoted (e.g., |BomlO[ p. 29]); this is the point to which 
the project had got by the time of Gourdon and Demichel's announcement. We- 
deniwski asserts in private communication that the project verified the first 10 12 
zeroes, and that the computation was double-checked (by the same method). 

The strongest claim prior to ZetaGrid was that of van de Lune (H = 3.293 TO 9 , 
first 10 10 zeroes; unpublished). Recently, Piatt [Pla| checked the first 1.1 • 10 11 
zeroes {H = 3.061 • 10 10 ) rigorously, following a method essentially based on that 
in [B00O6J. Note that [PlaJ uses interval arithmetic, which is highly desirable for 
floating-point computations. 

As for the verification of the binary Goldbach conjecture, we will use the work 
of Oliveira e Silva [OeS12j, which (as of April, 2012) states that every even integer 
4 < n < 4 • 10 18 is the sum of two primes. Presumably, this verification involves 
only integer arithmetic, and thus there are many fewer potential concerns than 
in the above. 
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Proposition A.l. Every odd integer 5 < n < uq is the sum of three primes, 
where 

'5.90698 • 10 29 if [GD04] is used - H = 2.44 • 10 12 ; 
n = < 6.15697 • 10 28 if ZetaGrid results are used (H = 2.419 • 10 n ), 
k 1.23163 • 10 27 if [Pli] is used ( H = 3.061 • 10 10 J. 

Proof. For n < 4-10 18 + 3, this is immediate from |OeS12j . Let 10 18 + 3 < n < no- 
We need to show that there is a prime p in [n — 4 — (n — 4)/A,n — 4], where 
A is large enough for (n — 4)/A < 4 • 10 18 — 4 to hold. We will then have that 
4 < n — p < 4 + (n — 4) / A < 4 • 10 18 . Since n — p is even, [QeS12j will then imply 
that n — p is the sum of two primes p', p" , and so 

n = p + p' + p" . 

Since n - 4 > 10 11 , the interval [n - 4 - (n - 4)/A, n - 4] with A = 28314000 
must contain a prime [RS03] . This gives the solution for (n — 4) < 1.1325 • 10 26 , 
since then (n - 4) < 4 • 10 18 - 4. Note 1.132 5 • 10 26 > e 59 . 

From here onwards, we use the tables in jHanul to find acceptable values of 
A. Since n — 4 > e 59 , we can choose 

'52211882224 if [GD04] is used (case (a)), 
A = < 13861486834 if ZetaGrid is used (case (b)), 
307779681 if (Pi] is used (case (c)). 

This gives us (n — 4) /A < 4 • 10 18 — 4 for n — 4 < e r °, where ro = 67 in case (a), 
ro = 66 in case (b) and ro = 62 in case (c). 

If n — 4 > e r ° , we can choose (again by [Ramd| ) 

'146869130682 in case (a), 
A = < 15392435100 in case (b), 
307908668 in case (c). 

This is enough for n — 4 < e 68 in case (a) , and without further conditions for (b) 
or (c). 

Finally, if n — 4 > e 68 and we are in case (a), [Ramd| assures us that the choice 

A = 147674531294 

is valid; we verify as well that (n - 4)/ A < 4 • 10 18 - 4. □ 

In the main theorem, I have chosen to include the assumption x > 3.1 • 10 28 . 
This means that, for it to be applied to Goldbach's three-prime problem (as in 
£j6]), we must assume uq is at least about 1.8 • 3.1 • 10 28 = 5.58 • 10 28 . In other 
words, we are using the ZetaGrid results (or something slightly weaker), but we 
never need [GD04] . 

Downgrading the assumptions so that only [Pla| is used is a fairly easy matter: 
we only need to modify some parameters starting in 15.2.11 The cost would be 
a moderate degradation in the bounds (at least about 5 percent) and a slight 
lengthening of the proof. 
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